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P : ABSTRACT 



We develop a new theoretical model for the spectral formation process in 
accretion-powered X-ray pulsars based on a detailed treatment of the bulk and 
c§ ! thermal Comptonization occurring in the accreting, shocked gas. A rigorous 

eigenfunction expansion method is employed to obtain the analytical solution for 
^ ! the Green's function describing the scattering of radiation injected into the col- 

umn from a monochromatic source located at an arbitrary height above the stellar 
surface. The emergent spectrum is calculated by convolving the Green's function 
with source terms corresponding to bremsstrahlung, cyclotron, and blackbody 
emission. The energization of the photons in the shock, combined with cyclotron 
absorption, naturally produces an X-ray spectrum with a relatively flat contin- 
uum shape and a high-energy quasi-exponential cutoff. We demonstrate that the 
new theory successfully reproduces the phase-averaged spectra of the bright pul- 
sars Her X-1, LMC X-4, and Gen X-3. In these luminous sources, it is shown that 
the emergent spectra are dominated by Gomptonized bremsstrahlung emission. 
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1. INTRODUCTION 

More than 100 pulsating X-ray sources have been discovered in the Galaxy and the Mag- 
ellanic Clouds since the seminal detection of Her X-1 and Cen X-3 more than three decades 
ago (Giacconi et al. 1971; Tananbaum et al. 1972). These sources display luminosities in the 
range Lx ~ io^^~^^ ergs s~^ and pulsation periods 0.1 s < P < 10^ s, and comprise a variety 
of objects powered by rotation or accretion, as well as several anomalous X-ray pulsars whose 
fundamental energy generation mechanism is currently unclear. The brightest pulsars are 
accretion-powered sources located in binary systems, with emission fueled by mass transfer 
from the "normal" companion onto one or both of the neutron star's magnetic poles. The 
accretion flow is channeled by the strong magnetic field into a columnar geometry, and the 
resulting emission is powered by the conversion of gravitational potential energy into kinetic 
energy, which escapes from the column in the form of X-rays as the gas decelerates through 
a radiative shock and settles onto the stellar surface. The spectra of the accretion-powered 
sources are often well fitted using a combination of a power-law spectrum in the 5 — 20 keV 
energy range plus a blackbody component with a temperature in the range T ~ 10^~^ K 
(e.g., Coburn et al. 2002; di Salvo et al. 1998; White et al. 1983) and a quasi-exponential 
cutoff at energy ii^ ~ 20 — 30 keV. There are also indications of cyclotron features and iron 
emission lines in a number of sources (e.g., Pottschmidt et al. 2005). The observations 
suggest typical magnetic field strengths ~ 10^^~^^ G. 

Previous attempts to calculate the spectra of accretion-powered X-ray pulsars based on 
static or dynamic theoretical models have generally yielded results that do not agree very well 
with the observed profiles (e.g., Meszaros & Nagel 1985a,b; Nagel 1981; Yahel 1980; Klein 
et al. 1996). Due to the lack of a fundamental physical model, most X-ray pulsar spectral 
data have traditionally been fitted using multicomponent functions of energy that include 
absorbed power-laws, cyclotron features, iron emission lines, blackbody components, and 
high-energy exponential cutoffs. The resulting model parameters are often difficult to relate 
to the physical properties of the source. However, the situation improved recently with the 
development by Becker & Wolff (2005a,b) of a new model for the spectral formation process 
based on the "bulk" or "dynamical" Comptonization (i.e., first-order Fermi energization) of 
photons due to collisions with the rapidly compressing gas in the accretion column . In the 
"pure" bulk Comptonization model of Becker & Wolff, the transfer of the proton kinetic 
energy to the photons, essential for producing the emergent spectrum, is rigorously modeled 
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using a transport equation that includes a first-order Fermi term that accounts for the energy 
gain experienced by photons as they scatter back and forth across the accretion shock. The 
seed photons for the upscattering process are provided by the blackbody radiation injected 
at the surface of the dense "thermal mound," located at the bottom of the accretion column 
(Davidson 1973). By ignoring the effect of thermal Comptonization in the column, the model 
corresponds physically to the accretion of a gas with thermal velocity much less than the 
dynamical (bulk) velocity. This "cold plasma" criterion is well satisfied in X-ray pulsars 
because the typical infiow speed at the top of the column is ~ 0.5 c and the characteristic 
temperature is ~ 10^ K. The resulting X-ray spectrum is characterized by a power law with 
photon index ctx > 2 combined with a low-energy turnover representing unscattered Planck 
radiation. The X-ray spectrum computed using the pure bulk Comptonization model has 
no high-energy cutoff, and therefore the photon index must exceed two in order to avoid an 
infinite photon energy density. Despite this restriction, the bulk Comptonization model has 
successfully reproduced the observed spectra of several steep-spectrum sources, including X 
Persei and GX 304-1 (Becker & Wolff 2005a,b). 

Many of the brightest X-ray pulsars, such as Her X-1 and Cen X-3, display spectra 
with photon indices ax < 2 in the 5 — 20keV energy range, combined with high-energy 
quasi-exponential cutoffs at 20 — 30keV (e.g.. White et al. 1983). This type of spectral 
shape cannot be explained using the pure bulk Comptonization model of Becker & Wolff 
(2005b). The presence of the exponential cutoffs at high energies, combined with the rela- 
tively flat shape at lower energies, suggests that thermal Comptonization is also playing an 
important role in luminous X-ray pulsars by transferring energy from high- to low-frequency 
photons via electron scattering. The two steps in the thermal process (Compton scattering 
followed by inverse-Conipton scattering) are described mathematically by the Kompaneets 
(1957) equation (see also Becker 2003). Although it is clear that the majority of the photon 
energization in X-ray pulsars occurs via the first-order Fermi (i.e., bulk Comptonization pro- 
cess), we find that in bright sources such as Her X-1, the role of the thermal Comptonization 
process must also be considered in order to reproduce the observed spectra. 

In this paper, we extend the pure bulk Comptonization model developed by Becker 
& Wolff (2005b) to include both bulk and thermal Comptonization by incorporating the 
full Kompaneets operator into the transport equation used to model the development of the 
emergent radiation spectrum. The exact solution to this equation yields the Green's function 
for the problem, which represents the contribution to the observed photon spectrum due to a 
monochromatic source at a fixed height in the accretion column. By exploiting the linearity 
of the mathematical problem, the Green's function can be used to compute the emergent 
spectrum due to any spatial-energetic distribution of photon sources in the column. We 
calculate X-ray pulsar spectra by convolving the Green's function with bremsstrahlung. 
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cyclotron, and blackbody sources, with the first two distributed throughout the column, and 
the latter located at the surface of the thermal mound. The accretion/emission geometry is 
illustrated schematically in Figure 1. Seed photons produced inside the column by various 
emission mechanisms experience electron scattering as they diffuse through the column, 
eventually escaping through the walls to form the emergent X-ray spectrum. The escaping 
photons carry away the kinetic energy of the gas, thereby allowing the plasma to settle onto 
the surface of the star. 

The remainder of the paper is organized as follows. In § 2 we briefly review the nature 
of the primary radiation transport mechanisms in the accretion column with a focus on dy- 
namical, thermal, and magnetic effects. The transport equation governing the formation of 
the radiation spectrum is introduced and analyzed in § 3, and in § 4 the exact analytical 
solution for the Green's function describing the radiation distribution inside the accretion 
column is derived. The spectrum of the radiation escaping through the walls of the accretion 
column is developed in § 5, and the physical constraints for the various model parameters 
are considered in § 6. The nature of the source terms describing the injection of blackbody, 
cyclotron, and bremsstrahlung seed photons into the accretion column is discussed in § 7. 
Emergent X-ray spectra are computed in § 8, and the results are compared with the ob- 
servational data for several luminous X-ray pulsars. The implications of our work for the 
production of X-ray spectra in accretion-powered pulsars are discussed in § 9. 



2. RADIATIVE PROCESSES 

The dynamics of gas accreting onto the magnetic polar caps of a neutron star was 
considered by Basko & Sunyaev (1976) and Becker (1998). The formation of the emergent 
X-ray spectrum in this situation was discussed by Becker & Wolff (2005b) based on the 
geometrical picture illustrated in Figure 1. Physically, the accretion scenario corresponds 
to the flow of a mixture of gas and radiation inside a magnetic "pipe" that is sealed with 
respect to the gas, but is transparent with respect to the radiation. The accretion column 
incorporates a radiation-dominated, radiative shock located above the stellar surface. Seed 
photons produced via a combination of cyclotron, bremsstrahlung, and blackbody radiation 
processes are scattered in energy due to collisions with electrons that arc infalling with high 
speed and also possess a large stochastic (thermal) velocity component. Blackbody seed 
photons are produced at the surface of the dense "thermal mound" located at the base of 
the flow, where local thermodynamic equilibrium prevails, and cyclotron and bremsstrahlung 
seed photons are produced in the optically thin region above the thermal mound. Hence the 
surface of the mound represents the "photosphere" for photon creation and absorption, and 
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the opacity is dominated by electron scattering above this point. 



2.1. Magnetic Effects 

The flow of gas in the accretion column of an X-ray pulsar is channeled by the strong 
{B ~ 10^^ G) magnetic fleld, and the presence of this field also has important consequences 
for the photons propagating through the plasma. In particular, vacuum polarization leads to 
birefringent behavior that gives rise to two linearly polarized normal modes (Ventura 1979; 
Nagel 1980; Chanan et al. 1979). The ordinary mode is polarized with the electric field 
vector located in the plane formed by the pulsar magnetic field and the photon propagation 
direction. For the extraordinary mode, the electric vector is oriented perpendicular to this 
plane. The nature of the photon-electron scattering process is quite different for the two po- 
larization modes, and it also depends on whether the photon energy, e, exceeds the cyclotron 
energy, e^, given by 

e, = ^ 11.57 Bu keV , (1) 

ZTrnieC 

where B12 = i?/(10^^G) and c, h, nie, and e represent the speed of light, Planck's con- 
stant, and the electron mass and charge, respectively. Ordinary mode photons interact with 
electrons via continuum scattering, governed by the approximate cross section (Arons et al. 
1987) 

c^ordle, V^) = <^ + k{e) cos^ ^p] , (2) 

where ip is the propagation angle with respect to the magnetic field, is the Thomson cross 
section, and 

^6) ^ I I ' 2 ' ^ ' (3) 

The scattering cross section for the extraordinary mode photons is more complex because 
these photons interact with the electrons via both continuum and resonant processes. In this 
case the total cross section can be approximated using (Arons et al. 1987) 

(7ext(e, <P) = cr^k{6) + ai (j)i{e, (fi) , (4) 

where 0; is the unity normalized line profile function (which is resonant at the cyclotron 
energy Ec), and 

(7; = 1.9 X 10^(7^ Sf2^ . (5) 

Equation (4) indicates that the extraordinary mode photons experience a substantial en- 
hancement in the scattering cross section close to the cyclotron energy, which refiects their 
ability to cause radiative excitation of electrons from the ground state to the first excited 
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Landau level (Ventura 1979). The excitation process is almost always followed by a radiative 
deexcitation, and therefore cyclotron absorption can be viewed as a form of resonant scat- 
tering (Nagel 1980). The resonant nature of the cyclotron interaction gives rise to a strong 
high-energy absorption feature in many observed X-ray pulsar spectra. 

The energy and angular dependences of the electron scattering cross section are quite 
different for the two polarization modes. For radiation propagating perpendicular to the 
magnetic field with energy e < ec, photons in the ordinary mode will see a cross section that is 
essentially Thomson, whereas the extraordinary mode photons will experience a cross section 
that is reduced from the Thomson value by the ratio (e/ec)^. On the other hand, in the limit 
of parallel or antiparallcl propagation relative to the field direction, both the ordinary and 
extraordinary mode cross sections are reduced by the factor (e/ec)^ relative to the Thomson 
value if e < Ec- In practice, the two modes communicate via mode conversion, which occurs 
at the continuum scattering rate (Arons et al. 1987). Hence if e is sufficiently below the 
cyclotron energy ec so that the resonant contribution to the cross section in equation (4) is 
negligible, then the scattering of photons propagating perpendicular to the magnetic field is 
dominated by the ordinary mode cross section, Uord = cTt; because this yields the smallest 
mean free path. 

Due to the importance of radiation pressure and the complexity of the electron scattering 
cross sections, the dynamical structure of the fiow is closely tied to the spatial and energetic 
distribution of the radiation. The coupled radiation-hydrodynamical problem is so complex 
and nonlinear that it is essentially intractable. In order to make progress a set of simplifying 
assumptions must be adopted. In particular, a detailed consideration of the angular and 
energy dependences of the electron scattering cross sections for the two polarization modes 
is beyond the scope of the present paper. We shall therefore follow Wang & Prank (1981) 
and Becker (1998) by treating the directional dependence of the electron scattering in an 
approximate way in terms of the constant, energy- and mode-averaged cross sections (7\\ and 
(T_L describing respectively the scattering of photons propagating either parallel or perpen- 
dicular to the magnetic field. For the magnetic field strengths B ~ 10^^^^'^ G and electron 
temperatures Tg ~ 10^^ ''K typical of luminous X-ray pulsars, the mean photon energy, e, is 
well below Cc- In this case the mean scattering cross section for photons propagating parallel 
to the field can be approximated by writing (see eqs. [2] and [4]) 



and the mean scattering cross section for photons propagating perpendicular to the field 
(which is dominated by the ordinary mode) is given by 




(6) 



(7_L ~ (7^ 



(7) 
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Note the substantial decrease in the opacity of the gas along the magnetic field indicated by 
equation (6), which is consistent with results presented by Canuto et al. (1971). 

In our calculations, the perpendicular scattering cross section is set equal to the 
Thomson value using equation (7). However, the value of the parallel scattering cross section 
(7|| is more problematic since the mean energy e cannot be calculated until the radiative 
transfer problem is solved. In practice, the value of g\\ is calculated using a dynamical 
relationship based on the flow of radiation-dominated gas in a pulsar accretion column, as 
discussed in § 6, where we also examine the self-consistency of our results. 



2.2. Radiation-Dominated Flow 

Radiation pressure governs the dynamical structure of the accretion fiows in bright 
pulsars when the X-ray luminosity Lx satisfies (Becker 1998; Basko & Sunyaev 1976) 



^crit 



' ' -— ergs s , (8) 



^/oTo^ \R 



where Tq is the polar cap radius, and denote the stellar mass and radius, respectively, 
and Lx is related to the accretion rate M via 

L. = '^. (9) 

When the luminosity of the system is comparable to Lcrit, the radiation fiux in the column 
is super-Eddington and therefore the radiation pressure greatly exceeds the gas pressure 
(Becker 1998). In this situation the gas passes through a radiation-dominated shock on its 
way to the stellar surface, and the kinetic energy of the particles is carried away by the high- 
energy radiation that escapes from the column. The strong gradient of the radiation pressure 
decelerates the material to rest at the surface of the star. The observation of many X-ray 
pulsars with Lx ~ 10'^^^^^ ergs s^^ implies the presence of radiation-dominated shocks close 
to the stellar surfaces in these systems (White et al. 1983; White et al. 1995). Note that 
radiation-dominated shocks are continuous velocity transitions, with an overall thickness of 
a few Thomson scattering lengths, unlike traditional (discontinuous) gas-dominated shocks 
(Blandford k, Payne 1981b). Despite the fact that the luminosities of the brightest X-ray 
pulsars exceed the Eddington limit for a neutron star, the material in the accretion column 
is decelerated to rest at the stellar surface, rather than being blown away, because the 
scattering cross section o"|| for photons propagating parallel to the magnetic field is generally 
much smaller than the Thomson value. 
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Becker & Wolff (2005b) considered tfie effects of pure bulk Comptonization on the 
emergent spectrum in an X-ray pulsar based on the exact velocity profile derived by Becker 
(1998) and Basko & Sunyaev (1976), which describes the accretion of gas onto a neutron star 
through a standing, radiation-dominated shock. In the present paper, our goal is to extend 
the model to include the effects of thermal Comptonization, which, along with cyclotron ab- 
sorption, is expected to produce a steepening of the spectrum at high energies, as observed 
in many luminous X-ray pulsars. The thermal process will also cause a flattening of the 
spectrum at lower energies, relative to the spectrum resulting from pure bulk Comptoniza- 
tion, due to the redistribution of energy from high- to medium-energy photons via electron 
recoil. Since there are many physical effects involved in the scenario considered here, this is 
clearly a rather complex problem. In order to render the problem mathematically tenable, 
we will utilize an approximate velocity profile that agrees reasonably well with the exact 
profile derived by Becker (1998) and used by Becker & Wolff (2005a,b) in their study of pure 
bulk Comptonization in X-ray pulsars. The validity of this approximation will be carefully 
examined by comparing the spectrum obtained using the approximate velocity profile with 
that computed using the pure bulk Comptonization model of Becker & Wolff in the limit of 
zero electron temperature, in which case both models should agree. 



2.3. Thermal Equilibration 

Another important question concerns the nature of the energy distribution in the accret- 
ing plasma. In an X-ray pulsar accretion column, the electrons will have an anisotropic energy 
distribution, described by a one-dimensional Maxwellian with a relatively high temperature 
(Tg ~ 10^ K) along the direction of the magnetic field, and an essentially monoenergetic 
distribution in the perpendicular direction (e.g., Arons et al. 1987). Most of the electrons 
reside in the lowest Landau state (the ground state), and these particles possess no gyrational 
motion. However, the electrons in the first excited Landau level have energy ec ~ 20 — 40keV 
(sec eq. [1]), which is much larger than the typical thermal energy in the parallel direction. 
The protons, which are not as strongly effected by the field, have a three-dimensional ther- 
mal distribution. Collisions with high speed protons can cause the electrons to be excited 
to higher Landau levels, followed rapidly by radiative deexcitation. This process essentially 
converts proton kinetic energy into radiation energy via the production of cyclotron photons, 
which cools the plasma. The Comptonization of these seed photons leads to further cooling 
before the radiation escapes through the walls of the accretion column in the form of X-rays. 

In principle, the temperature of the ions may depart significantly from the electron tem- 
perature, depending on the ratio of the Coulomb coupling rate to the dynamical timescale 



-9- 



for accretion onto the stellar surface. We can determine if the gas in the X-ray pulsar accre- 
tion column has a two-temperature structure by comparing the electron-proton equilibration 
timescale with the dynamical and radiative cooling timescales. Based on equation (103) from 
Arons et al. (1987), we can approximate the electron-ion coupling timescale using 

" rp \ 3/2 



where p is the density of the accreting gas. In the bright sources of interest here, we typically 
find that Tg ~ 10^"^ K and p ~ 10~^ — 10~^gcm~^ within the radiating portion of the 
accretion column, giving t^i ~ 10~^ — 10~^^ s. 

The opacity of the gas in the accretion column of a bright X-ray pulsar is dominated 
by electron scattering, and therefore inverse- Compton scattering of soft photons by the hot 
electrons is the primary cooling mechanism. The rate of change of the photon energy density 
due to inverse- Compton scattering is given by (see eq. [7.22] from Rybicki & Lightman 1979) 



dU 
'dt 



- ^^^^ TT (^^^ 
riecrc U , (11) 



IC TTlgC 



where Ug = p/rrip is the electron number density for pure, fully-ionized hydrogen, nip is 
the proton mass, U denotes the photon energy density, and a represents the angle-averaged 
electron scattering cross section, which depends on the energy and angular distribution of 
the radiation field. In general, we expect to find that a\\ < a < a±, where a\\ and a± are 
the scattering cross sections for photons propagating either parallel or perpendicular to the 
magnetic field (see eqs. [6] and [7]). 

The photon energy density U can be estimated in terms of the X-ray luminosity Lx and 
the column radius Tq using 

Ur^J^. (12) 

TTTq c 

By combining equation (11) and (12), we find that the inverse-Compton cooling timescale 
for the electrons is given by 

^^^= dU/dt [T^J yik^J Uo^^ergss-iJ ^^^^ 

In the X-ray pulsar application, we typically obtain (see § 6) o'/a^ ~ 10~^ — 10~^, Tq ~ 0.1 — 
1 km, and Lx ^ 10^^ ergs s~^, and therefore tic ~ 10~^ s. Next we recall that the characteristic 
dynamical (free-fall) timescale for neutron star accretion is given by 
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thermal mound polar cap 

Fig. 1. — Schematic depiction of gas accreting onto the magnetic polar cap of a neutron star. 
Seed photons are created throughout the column via bremsstrahlung and cyclotron emission, 
and additional blackbody seed photons are emitted from the surface of the thermal mound, 
near the base of the column. 

It is apparent that the dynamical timescale tdyn and the inverse-Compton cooling timescale 
tic each exceed the electron-ion equilibration timescale tei by several orders of magnitude, 
and therefore we conclude that the electron and ion temperatures are essentially equal in 
the sources of interest here. 

2.4. Thermal and Bulk Comptonization 

The strong compression that occurs as the plasma crosses the radiative shock renders 
it an ideal site for the bulk Comptonization of seed photons produced via bremsstrahlung, 
cyclotron, or blackbody emission. In the bulk Comptonization process, particles experience 
a mean energy gain if the scattering centers they collide with are involved in a converging 
flow (e.g., Laurent & Titarchuk 1999; TuroUa, Zane, & Titarchuk 2002). By contrast, in the 
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thermal Comptonizaton process, particles gain energy due to the stochastic motions of the 
scattering centers via the second-order Fermi mechanism (e.g., Sunyaev & Titarchuk 1980; 
Becker 2003). In the X-ray pulsar application, the scattering centers are infalling electrons, 
and the energized "particles" are photons. Since the inflow speed of the electrons in an X- 
ray pulsar accretion column is much larger than their thermal velocity, bulk Comptonization 
dominates over the stochastic process except at the highest photon energies (Titarchuk, 
Mastichiadis, & Kylafis 1996). 

The fundamental character of the Green's function describing both thermal and bulk 
Comptonization was studied by Titarchuk & Zannias (1998) for the case of accretion onto a 
black hole. These authors established that the Green's function can be approximated using 
a broken power-law form with a central peak between the high- and low-energy portions of 
the spectrum, for either type of Comptonization. Furthermore, they concluded that bulk 
Comptonization dominates over the thermal process if the electron temperature Te < 10^ K. 
The theory developed in the present paper represents an extension of the same idea to neu- 
tron star accretion, resulting in the exact Green's function for the X-ray pulsar spectral 
formation process. In agreement with Titarchuk & Zannias (1998), we find that the bulk 
Comptonization process dominates the energy exchange between the electrons and the pho- 
tons for the temperature range relevant in X-ray pulsars. However, we also establish that 
thermal Comptonization nonetheless plays a central role in forming the characteristic spec- 
tral shapes noted in the bright sources. For example, Compton scattering of the high-energy 
photons is crucial for producing the observed quasi-exponential turnovers at high energies as 
a result of electron recoil, which is part of the thermal Comptonization process represented 
by the Kompaneets operator discussed in § 3. The subsequent inverse-Compton scattering of 
soft photons by the recoiling electrons redistributes the energy to lower frequency photons, 
helping to flatten the spectra as observed. 

Lyubarskii & Sunyaev (1982) analyzed the effect of bulk and thermal Comptonization 
in the context of neutron star accretion using the same dynamical model for the flow con- 
sidered here. However, they did not include the effect of photon escape from the accretion 
column, which renders their model inapplicable for the modeling of spectral formation in 
X-ray pulsars. In this paper we extend the work of Lyubarskii & Sunyaev by including the 
crucial effect of photon escape. Poutanen & Gierhnski (2003) computed X-ray pulsar spectra 
based on the thermal Comptonization of soft radiation in a hot layer above the magnetic 
pole, but their model did not include a complete treatment of the bulk process, which is 
of central importance in X-ray pulsars. The new theory presented here therefore represents 
the flrst exact, quantitative analysis of the role of bulk and thermal Comptonization in the 
X-ray pulsar spectral formation process. 



3. FORMATION OF THE RADIATION SPECTRUM 



We follow the approach of Basko & Sunyaev (1975, 1976) and Becker & Wolff (2005b), 
and assume that the upstream flow is composed of pure, fully-ionized hydrogen moving 
at a highly supersonic speed, which is the standard scenario for accretion-powered X-ray 
pulsars. Our transport model employs a cylindrical, plane-parallel geometry, and therefore 
the velocity, density, and pressure are functions of the distance above the stellar surface, but 
they are all constant across the column at a given height (see Fig. 1). In the region above 
the thermal mound in a luminous X-ray pulsar, the gas is radiation-dominated, and the 
photons interact with the matter primarily via electron scattering, which controls both the 
spatial transport and the energization of the radiation (Arons et al. 1987). "Seed" photons 
injected into the flow are unable to diffuse very far up into the accreting gas due to the 
extremely high speed of the inflow. Most of the photons therefore escape through the walls 
of the column within a few scattering lengths of the mound, forming a "fan" type beam 
pattern, as expected for accretion-powered X-ray pulsars (e.g., Harding 1994, 2003). We 
are interested in obtaining the steady-state, polarization mode-averaged photon distribution 
function, f{z,e), measured at altitude z and energy e inside the column resulting from the 
reprocessing of blackbody, cyclotron, and bremsstrahlung seed photons. The normalization 
of / is deflned so that e^f{z, e) de gives the number density of photons in the energy range 
between e and e -|- de, and therefore / is related to the occupation number distribution n via 



In the cylindrical, plane-parallel geometry employed here, the photon distribution f{z, e) 
satisfles the transport equation (e.g., Becker & Begelman 1986; Blandford & Payne 1981a; 



where z is the distance from the stellar surface along the column axis, v < is the inflow 
velocity, Q denotes the photon source distribution, and t^sc represents the mean time photons 
spend in the plasma before diffusing through the walls of the column. The source function Q 
is normalized so that Q{z, e) de dz gives the number of seed photons injected per unit time 
between z and z + dz with energy between e and e -|- de. The left-hand side of equation (15) 



/ = 87rn/(c3/i3). 



3.1. Transport Equation 



Becker 2003) 




(15) 
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denotes the comoving time derivative of the radiation distribution /, and the terms on the 
right-hand side represent first-order Fermi energization ("bulk Comptonization"), spatial 
diffusion along the column axis, photon escape, thermal Comptonization, and photon injec- 
tion, respectively. We are interested here in the steady-state version of equation (15) with 
df/dt = 0. The transport equation employed here is similar to the one analyzed by Becker & 
Wolff (2005b), except that thermal Comptonization has now been included via the appear- 
ance of the Kompaneets (1957) operator, and the source term Q has been generalized to treat 
the production of seed radiation throughout the column, rather than focusing exclusively on 
the injection of blackbody photons at the surface of the thermal mound. Equation (15) does 
not include the term used by Becker & Wolff (2005b) to describe absorption at the thermal 
mound surface because this effect is negligible when bremsstrahlung and cyclotron emission 
are included in the model, as discussed in § 8. The total photon number and energy densities 
associated with the radiation distribution / are given, respectively, by 

/•OO /"OO 

n{z)^ e^f{z,e)de, U{z) = e^f{z,e)de. (16) 
JO Jo 

Following Becker (1998) and Becker & Wolff (2005b), we compute the mean escape time 
using the diffusive prescription 

^esc(^) = , r±{z) = nea±ro , (17) 

c 

where t± represents the perpendicular scattering optical thickness of the cylindrical accretion 
column, and t± and ^esc a-re each functions of z through their dependence on the electron 
number density rig. Becker (1998) confirmed that the diffusion approximation employed in 
equation (17) is valid since t± > 1 for typical X-ray pulsar parameters. Equation (17) for 
the mean escape timescale can be rewritten as 

tesciz) = ^ , 18 

Tririp c \v\ 

where 

M = nrl p \v\ — constant , (19) 

and p — nip Ug. Since the escape timescale is inversely proportional to the flow velocity, the 
column becomes completely opaque at the surface of the neutron star due to the divergence 
of the electron number density there. The relationship between the escape-probability for- 
malism employed here and the physical distribution of radiation inside the accretion column 
was discussed in detail by Becker & Wolff (2005b) in their § 6.3. 

In our approach to solving equation (15) for the spectrum f{z, e) inside the accretion 
column, we shall first obtain the Green's function, f^{zo, z,eo,e), which is the radiation 
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distribution at location z and energy e resulting from the injection of TVq photons per second 
with energy eo from a monochromatic source at location zq. The determination of the Green's 
function is a useful intermediate step in the process because it provides us with fundamental 
physical insight into the spectral redistribution process, and it also allows us to calculate the 
particular solution for the spectrum / associated with an arbitrary photon source Q{z^e) 
using the integral convolution (Becker 2003) 

/(z, e) = r r el Q{zo, 60) deo dzo . (20) 

^0 ^0 Wo 

The technical approach used to solve for the Green's function, carried out in § 4, involves the 
derivation of eigenvalues and associated eigenfunctions based on the set of spatial boundary 
conditions for the problem (see, e.g., Blandford & Payne 1981b; Payne & Blandford 1981; 
Schneider & Kirk 1987; Colpi 1988). 

The steady-state transport equation governing the Green's function is (cf. eq. [15]) 



G 



dz dz 3 de dz \3nea|| dz J t. 

n^c 1 d 



+ 



rrieC^ e^ de 



^ NoS{e-eo)6{z- zq) ^^i) 



de 

and the associated radiation number and energy densities are given by (cf. eqs. [16]) 

/•OO /"OO 

^g(^) = / e'^ fQ{zQ,z,eQ,e)de , U^{z)= / e^ f^{zQ,z,eQ,e) de . (22) 

JQ Jo 

Following Lyubarskii & Sunyaev (1982), we will assume that the electron temperature has 
a constant value, which is physically reasonable since most of the Comptonization occurs in 
a relatively compact region above the thermal mound. However, it is important to confirm 
the validity of this assumpton using a detailed numerical model that incorporates a self- 
consistent calculation of the temperature distribution, which wc plan to develop in future 
work. Under the assumption of a constant electron temperature, it is convenient to work in 
terms of the dimensionless energy variable x, defined by 

X(e) . ^ . (23) 

We can make further progress by transforming the spatial variable from z to the scattering 
optical depth parallel to the magnetic field, r, which is related to z via 



dr 



n(.{z)a\\dz, t{z) = / ne{z')a\idz' , (24) 

Jo 
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so that z and r both vanish at the stellar surface. 

Making the change of variable from z to r in equation (21), we find after some algebra 
that the transport equation for the Green's function can be written in the form 

vdj^ _ Idvxdf l^_ejlr 



c dr c dr 3 dx 3 ° 

+ 



a kTe 1 d 



an rrieC^ x^ dx 



where Tq = t{zo), Xo = xi^o)^ we have introduced the dimensionless parameter 

^^nr^rn^ (26) 

which determines the importance of the escape of photons from the accretion column. Becker 
(1998) derived the exact solution for the flow velocity profile in a radiation-dominated pulsar 
accretion column, and demonstrated that the condition ^ = must be satisfied in order 

to ensure that the flow comes to rest at the stellar surface. This condition represents a 
balance between the characteristic timescale for photon escape and the dynamical (accretion) 
timescale, as discussed in § 6. Physically, this balance reflects the requirement that the 
kinetic energy of the flow must be radiated away through the column walls in the same 
amount of time required for the gas to settle onto the star. Note that we can write the 
Green's function as either /^{zq, z, cq, e) or /g(to, t, Xo, x) since the variables {z, e) and (x, r) 
are interchangeable via equations (23) and (24). 



3.2. Separability 

Lyubarskii & Sunyaev (1982) demonstrated that when x 7^ Xo, the transport equa- 
tion (25) is separable in energy and space if the velocity profile has the particular form 

v{t) — —a CT , (27) 

where a is a positive constant. By combining equations (19), (24), (26), and (27), we can 
express r as an explicit function of the altitude z, obtaining 

\ 1/4 / o \ 1/2 



Using this result to substitute for r in equation (27), we note that the velocity profile required 
for separability is related to z via 

\ 1/4 /o \ 1/2 

ail \ ' f 2az^ ' 
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Although this profile describes a fiow that stagnates at the stellar surface {t — 0, z — 0) 
as required, the details of the velocity variation deviate somewhat from the exact solution 
for the velocity profile in a radiation-dominated pulsar accretion column derived by Becker 
(1998) and Basko & Sunyaev (1976), which can be stated in terms of the altitude z using 



[z) = -V{{ 



1 - 



3 



(30) 



where the free-fall velocity from infinity onto the stellar surface, vg, and the altitude at the 
sonic point, ^sp; are given by 



-n^J . ..p.-s^-j in^.j. (31) 

In order to make further mathematical progress in the computation of the radiation 
spectrum emitted by the accretion column, we must utilize the approximate (separable) 
velocity profile given by equation (29). However, before doing so, we must ensure that the 
approximate velocity profile agrees reasonably well with the exact profile (eq. [30]) in the 
lower portion of the accretion column, where most of the spectral formation occurs. We can 
accomplish this by setting the two profiles equal to each other at the sonic point {z = Zgp), 
and using this condition to solve for the constant a, which yields 

49 1n(7/3) i?,c2 ' ^ > 

or, equivalently, 

a = ,.24^\(^\\. (33) 

This relation allows us to compute a as a function of ^ for given values of the stellar mass M,, 
and radius i?*. In radiation-dominated pulsar accretion columns, ^ = 2/\/3, and therefore 
a is of order unity (Becker 1998). 

The detailed shapes of the approximate (cq. [29]) and exact (cq. [30]) velocity profiles 
as functions of z arc compared in Figure 2. In general, the two functions agree fairly well, 
although the approximate profile overestimates the correct velocity close to the stellar sur- 
face. Nevertheless, since the overall compression associated with the approximate profile is 
correct, we expect that the net effect of bulk Comptonization will be accurately modeled 
using the approximate profile. We will therefore adopt the form for the velocity profile given 
in terms of r by equation (27) or in terms of z by equation (29), with a computed using 
equation (33). The validity of this approach will be carefully examined in § 6 by comparing 
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Fig. 2. — Comparison of the approximate {solid line) and exact {dashed line) velocity profiles 
evaluated using eqs. (29) and (30), respectively, in units of c. The constant a has been 
computed using eq. (32), which ensures that the two velocities agree at the sonic point, 

the results obtained for the Comptonization of radiation in a "cold" accretion flow (Tg 0) 
described by the approximate velocity profile with those computed using the model of Becker 
& Wolff (2005b), who treated the Comptonization of radiation in a cold accretion column, 
described by the exact velocity profile. 



4. EXACT SOLUTION FOR THE GREEN'S FUNCTION 

Adopting the velocity profile given by equation (27), we find that the transport equa- 
tion (25) for the Green's function can be reorganized to obtain 



3 dx mec2 dx 



XM /o + 



dfa 
dx 



1^ 
3 dr^ 



dr 
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+ 



nrlckTeel 



(34) 



Lyubarskii & Sunyaev (1982) analyzed this equation for the special case ^ = 0, which 
corresponds to the neglect of the escape of photons from the accretion column. Their results 
do not describe the formation of the emergent spectrum in an accretion-powered X-ray 

pulsar since photon escape is a critical part of that process. We therefore extend the results 
of Lyubarskii & Sunyaev below by deriving the closed-form solution for the Green's function 
in the X-ray pulsar problem for the general case with ^ > 0. 

When X Xo, the 5-function in the transport equation (34) makes no contribution, and 
therefore the differential equation is linear and homogeneous. The transport equation can 
then be separated in energy and space using the functions 



/a(t,x) = giKr) h{X,x) 



(35) 



where A is the separation constant. We find that the spatial and energy functions, g and h, 
respectively, satisfy the differential equations 



1 (fg dg 
3d^ + ""d^ + 



and 



1 d 

X^ dx 

where the parameter 8 is defined by 



dh 
dx 



-eVr2U = 0, 



-5x^-5\h^Q, 
dx 



(36) 



(37) 



^ _ Q; f7|| TOgC 



3 a kT, 



(38) 



This is equivalent to the quantity 5 introduced by Lyubarskii & Sunyaev (1982) if we set 
a — (7|| , since these authors did not include any angle dependence in the electron scattering 
cross section. 



4.1. Eigenvalues and Spatial Eigenfunctions 

In order to obtain the solution for the spatial separation function g, we must first 
consider the boundary conditions that g must satisfy. In the downstream region, as the 
gas approaches the stellar surface, we expect the advective and diffusive components of the 
radiation flux to vanish due to the divergence of the electron density. The advective flux 
is indeed negligible at the stellar surface since — > as r ^ (see eq. [27]). However, in 
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order to ensure that the diffusive fiux vanishes, we must require that dgi/dr— >Oasr— >0. 
Conversely, in the upstream region, we expect that g ^ as t ^ oo since no photons can 
diffuse to large distances in the direction opposing the plasma flow. With these boundary 
conditions taken into consideration, we find that the fundamental solution for the spatial 
separation function g has the general form 

f e-"(3+-)-'/^ M f a, i, ^) , T < To , 
9{\r)ocl , > J (39) 

where M and U denote confluent hypergeometric functions (Abramowitz & Stegun 1970), 
and we have made the definitions 

a.^!i±|^, (40) 

Note that in the radiation-dominated case with ^ 2/\/3, we obtain w — 5. 

Equation (36) is linear, second-order, and homogeneous, and consequently both the 
function g and its derivative dg/dr must be continuous at the source location, r = tq. 
The smooth merger of the M and U functions at the source location requires that their 
Wronskian, 

, , , / 1 awT'^\ d 1 awT^ 

u){t) = M a, -,—— —[/ a, - 



2' 2 ) dr \ T 2 
( 1 ciWT^\ d 1 awT^\ ,,,, 

- ^("•2— )5;^n°'2— ) ■ 

must vanish at r = Tq. This condition can be used to solve for the eigenvalues of the 
separation constant A. By employing equation (13.1.22) from Abramowitz & Stegun (1970) 
to evaluate the Wronskian, we obtain the eigenvalue equation 

.M = -£(lM(a^,.».V..O. (42) 
1 \a) 

The left-hand side vanishes when r(a) — > ±oo, which implies that a — —n, where n — 
0, 1, 2,. . . By combining this result with equation (40), we conclude that the eigenvalues A„ 
are given by 

Anw -|- w -|- 3 „ , , „, 

An = 2 ' ^ = 0'1'2,... (43) 

When A = A„, the spatial separation functions g reduce to a set of global eigenfunctions 
that satisfy the boundary conditions at large and small values of r. In this case, we can 
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use equations (13.6.9) and (13.6.27) from Abramowitz & Stegun to show that the confluent 
hypergcometric functions M and U appearing in equation (41) are proportional to the gen- 
eralized Laguerre polynomials and consequently the spatial eigenfunctions can be 
written as 

5„(r) ^ ^(A„, r) = e-"(3+-)^V4 4-1/2) (^^^ _ (44^ 

Based on equation (7.414.3) from Gradshteyn & Ryzhik (1980), we note that the spatial 
eigenfunctions satisfy the orthogonality relation 

foo r r(ra+l/2) _ 

/ e^-'/2^„(T) grnir) dr = , n -m , ^^^^ 

Jo 10, n f^m . 



4.2. Energy Eigenfunctions 

The solution for the energy separation function h depends on the boundary conditions 
imposed in the energy space. As e — > 0, we require that h not increase faster than 
since the Green's function must possess a finite total photon number density (see eq. [22]). 
Conversely, as e — > 00, we require that h decrease more rapidly than in order to ensure 
that the Green's function contains a finite total photon energy density. Furthermore, in 
order to avoid an infinite diffusive fiux in the energy space at x = xo, the function h must 
be continuous there. The fundamental solution for the energy eigenfunction that satisfies 
the various boundary and continuity conditions can be written as 

where M^^^ and W^^^ denote Whittaker functions, and we have made the definitions 

/.^ 1(5 + 4), /.^i[(3-(5)^ + 4a„]^/' . (47) 

Note that each of the eigenvalues A„ results in a different value for /x, and the parameter 5 
is defined in equation (38). Equation (46) can also be written in the more compact form 

hn{x) = X""'' e-^'^ M,,^(x^i„) W,,,{xra^) , (48) 



where 



Xmin = niin(x, Xo) , Xmax = max(x, Xo) 



(49) 
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4.3. Eigenfunction Expansion 

The spatial eigenfunctions Qniv) form an orthogonal set, as expected since this is a 
standard Sturm-Liouville problem. The solution for the Green's function can therefore be 
expressed as the infinite series 

oo 

/g(^o, r,Xo,X) = ^ Cngn{r) K{x) , (50) 

n=0 

where the expansion coefficients C„ are computed by employing the orthogonality of the 
eigenfunctions, along with the derivative jump condition 



dfo 



lim 

£^0 dx 



X=XO+e 



dx 



2774 ' (51) 



^=^,-e a Tiro ceo 



which is obtained by integrating the transport equation (34) with respect to x in a very 
small range surrounding the injection energy xo- Substituting using the expansion for 
yields 

lim E C„ gM [KiXo + e)- K{xo - e)] ^ - '''''^ ^Ifj' , (52) 

£^0 ^ — ' a ITTn C en 

where primes denote differentiation with respect to x- By employing equation (48) for hn, 
we find that 

CngniVXo e W{Xq) = 2 4 ' ^3) 

where the Wronskian, 2U, is defined by 

2n(xo) = WUxo) - W,,,(xo) K^(Xo) . (54) 

The Wronskian can be evaluated analytically using equations (13.1.22), (13.1.32), and 
(13.1.33) from Abramowitz & Stegun (1970), which yields 

Using this result to substitute for W{xo) in equation (53) and reorganizing the terms, we 
obtain 

r(l + 2//)C,^„(r) ^ 3NoSe^^/H{T-To) 
^ Fit,- K + 1/2) anrlcx^im' ' ^ ^ 

where is a function of A„ via equation (47). We can now calculate the expansion co- 
efficients C„ by utilizing the orthogonality of the spatial eigenfunctions gn represented by 
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equation (45). Multiplying both sides of equation (56) by e^"'^^/^g'^(r) and integrating with 
respect to r from zero to infinity yields, after some algebra, 

^ 3No 6 e^°/^ e^""° r{n-K + 1/2) n! gM 

Tirlc{kT,)\^,^ r(l + 2^)r(n + l/2) ' ^ > 

The final closed-form solution for the Green's function, obtained by combining equations (48), 
(50), and (57), is given by 

3A(,(W-^^'^-/^v^ ^'^-i,.(Yo-v)/2 °^ r(//-/v + l/2)//! 
/c(^o,r,Xo,Xj - ^rlc{kT,fx^V^ ^ r(l + 2/.) r(n + 1/2) 

X 9n{To) gn{r) M«,^(Xmin) W^K,M(Xmax) , (58) 

where the spatial eigenfunctions Qn are computed using equation (44) and the parameters k, 
and n are given by equations (47). The Green's function can also be expressed directly in 
terms of the photon energy e by writing 



3No5 kTe e3"^o/2^ ^.-4 ^{eo-e)/i2kT,) - r(/x - K + 1/2) n! 
/g('^o,t, eo,e) = 2 — ^-1= 



^rlcel^ ^ r(l + 2/x)r(n + l/2) 

X 9n{T,) 9M M.,, (g) W^,, , (59) 

where 

Cmin = min(e, eo) , e^ax = max(e, eo) . (60) 

This exact, analytical solution for is one of the main results of the paper, and it provides 
a very efficient means for computing the steady-state Green's function resulting from the 
continual injection of monochromatic seed photons from a source at an arbitrary location 
inside the accretion column. The eigcnfunction expansion converges rapidly and therefore 
we can generally obtain an accuracy of at least four significant figures in our calculations of 
/q by terminating the series in equations (58) or (59) after the first 5-10 terms. 



4.4. Asymptotic Povi^er-Lavi^ Behavior 

Asymptotic analysis of the energy function W^^/^ appearing in equation (59) reveals that 
if the photon energy e and injection energy eo satisfy the conditions eo < e <S kTg, then the 
spectrum has the power-law form 
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where 5 is given by equation (38) and Aq is the leading eigenvalue computed using equa- 
tion (43) with n = 0. We will demonstrate in § 6 that bulk Comptonization dominates over 
thermal Comptonization in the limit of large S, in which case equation (61) reduces to 

(X e-^° , 5^oo, (62) 

which is the same solution obtained by Becker & WolfT (2005b) in the case of pure bulk 
Comptonization. However, in contrast to the pure-bulk model, in the scenario considered 
here the power-law shape only extends up to the Wien turnover at photon energy e ~ kT^. 

We can gain some insight into the role of thermal Comptonization by comparing the 
power-law indices in equations (61) and (62). In general, we find that 

Ao > 3-^W(3-^)- + 4V _ (,3) 

with the equality holding in the hmit of large 5. This implies that for small values of 5, 
thermal Comptonization causes a flattening of the spectrum in the region e <^ kT^ due to 
the transfer of energy from high- to medium-energy photons via electron recoil. If the flow 
is radiation-dominated, as expected in the bright pulsars, then ^ = 2/-\/3 and we obtain 

w — h and Ao = 4 according to equations (40) and (61), respectively. Luminous pulsars 
are generally dominated by bulk Comptonization, and therefore it follows that the photon 
spectral index ax = Aq — 2 = 2 in the region of the spectrum below the Wien cutoff. 



5. SPECTRUM OF THE ESCAPING RADIATION 

Equations (58) and (59) represent the exact solution for the Green's function de- 
scribing the radiation spectrum inside a pulsar accretion column resulting from the injection 
of A'^o seed photons per unit time from a monochromatic source located at r = Tq (or, equiv- 
alently, at z = zq). Since the fundamental transport equation (15) is linear, wc can use 
the analytical solutions for the Green's function to calculate the spectrum of the radiation 
escaping from the accretion column for any desired source distribution. 



5.1. Green's Function for the Escaping Radiation Spectrum 

In the escape-probability approach employed here, the associated Green's function for 
the number spectrum of the photons escaping through the walls of the cylindrical column is 
computed using 

N^(zo, z, eo, e) = ° , /^(^o, z, eo, e) , (64) 
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where the escape timescale t^sc is evaluated as a function of z by combining equations (18) 
and (29), which yields 

te.o{z) = ■ 65 

\ zTT nip a\\az j 

The quantity dz de represents the number of photons emitted from the disk-shaped vol- 
ume between positions z and z -\- dz per unit time with energy between e and e -|- de. 

Since the quantities (^, ^o) and (r, tq) are interchangeable via equations (24), we are free 
to work in terms of the more convenient parameters (r, Tq) without loss of generality. In this 
case equations (26), (28), and (65) can be combined to reexpress the escape timescale as 

^esc(r) = ^ f^V' . (66) 



By using this result to substitute for ^esc in equation (64), we find that the Green's function 
for the escaping radiation spectrum is given in terms of the optical depth by 

\ 1/2 

(Til ^ ' 



A^, (To,T,eo,e) = Trroca^ 1^— J re /^(ro, r, eo, e) , (67) 

where is computed using equation (59). The factor of r on the right-hand side of equa- 
tion (67) indicates that the emitted radiation is strongly attenuated near the stellar surface 
due to the divergence of the electron number density, which inhibits the escape of the pho- 
tons through the walls of the accretion column. Equations (59) and (67) can be combined 
to express the closed-form solution for Nf as 

. „ 37Vo (5 C kT^ J2^i^\ e3"-<?/2e'^-2 f,{eo-e)/{2kn) ^ 
N^{to, t, eo, e) = ^ 

_ r(i + 2^)r(n + i/2) [W, ) [Wj' 



X 

n=0 



where e^ain and Cmax are defined by equations (60), and w, and ii are computed using 
equations (40) and (47). 

In Figure 3 we plot the Green's function for the escaping photon number distribution, 
as a function of the photon energy e and the scattering optical depth r above the 
stellar surface for two values of the parameters a and cq using equation (68). In order 
to clearly illustrate the general features of the spectrum, we set a\\ = lO^^cr^, (t±_ — cr^, 
a = 0.1 (7^, TVo = 1, C = 2/v/3, Tq = 0.1, Tq = 10^ cm, and = 10^ K. We show in § 6 that 
the parameter 5 defined by equation (38) determines the relative importance of bulk and 
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Fig. 3. — Green's function (s~^ cm~^ keV~^) describing the photon number spectrum 
escaping from an X-ray pulsar accretion column as a function of the radiation energy e 
(keV) and the scattering optical depth r above the surface, evaluated using eq. (68). The 
curves were computed by setting (T|| = 10"'^ a^, a± = a^, a = 0.1 a^, Nq = 1, ^ = 2/^/3, 
tq = 0.1, tq = 10"^ cm, and Tg = 10^ K. The values of r, eo, and a are indicated for each plot, 
and the values of 6 are discussed in the text. 

thermal Comptonization. The values of 6 obtained in Figure 3 are 6 = 0.20 for a = 0.1, 
and 5 = 0.79 for a = 0.4. For the cases with eo = 0.1 keV, a distinct Wien hump is visible 
at e ~ 1 keV due to thermal Comptonization. However, the hump is less pronounced when 
a = 0.4 due to the greater strength of bulk Comptonization, which tends to create a power- 
law spectral shape with photon index ax = 2 (see § 4.4). The Wien feature is almost invisible 
when Co = 1 keV because in this case the energy of the injected photons is comparable to 
the thermal energy of the electrons. Note that for small values of r, the escape of radiation 
is inhibited because of the divergence of the electron density near the stellar surface. Little 
radiation is able to escape from the column for r S> 1 because the photons are advectively 
trapped near the stellar surface due to the high-speed inflow (see § 6). Hence the spatial 
distribution of the escaping radiation peaks around r ~ 1. 
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5.2. Altitude-Dependent Spectrum for an Arbitreiry Source 

Since the fundamental transport equation (15) is linear, the analytical results for the 
Green's function obtained in § 5.1 provide the basis for the consideration of any source dis- 
tribution. By analogy with equation (64), we can compute the photon spectrum emitted 
through the walls of the accretion column for an arbitrary source distribution Q in equa- 
tion (15) by writing 

N^{z,e)^^^fiz,e) , (69) 

'■esc V '2' J 

where / is the particular solution calculated using the integral convolution given by equa- 
tion (20). Substituting for and / in equation (20) using equations (64) and (69), respec- 
tively, we find that the particular solution for the emitted photon spectrum can be written 

as 

K{z, e) = r r ^e{zo,z,eo,e) ^^^^^ _ ^^^^ 

Jo Jo Wq 

The spatial integration converges despite the infinite upper bound for zq because the seed 
photon distribution Q is locahzed in the lower region of the accretion column (see § 6.1). 
Equation (70) facilitates the calculation of the photon spectrum emitted by the accretion 
column as a function of photon energy e and altitude z for any source distribution Q. How- 
ever, due to the large distances to the known pulsars, current observations are unable to 
resolve the spatial distribution of the emission from the column. It is therefore necessary to 
integrate the emitted radiation field with respect to z in order to compare the theoretical 
predictions with the available spectral data, as we discuss below. 



5.3. Column-Integrated Escaping Green's Function 

By integrating over the vertical structure of the accretion column, we can compute 
the total emitted radiation distribution, which corresponds approximately to the phase- 
averaged spectrum of the X-ray pulsar. In the cylindrical geometry employed here, the 
formal integration domain is the region < 2; < z^a^, where z^a^ is the altitude at the 
upper surface of the radiating region within the accretion column, determined by setting the 
inflow velocity equal to the local free-fall velocity (see § 6.1). In practice, however, we can 
replace the upper integration bound ^^max with infinity without introducing any significant 
error because very few photons escape more than one or two scattering optical depths above 
the stellar surface (see § 6.5 and Becker 1998). The replacement of the upper bound is 
advantageous because the resulting integral can be performed analytically. For the case of 
a monochromatic source, we therefore define the column-integrated Green's function for the 
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escaping photon spectrum using 

POD 

^f{zo,eo,e)= / N^{zo,z,eo,e)dz , (71) 
Jo 

where $f de represents the total number of photons escaping from the column per unit time 
with energy between e and e + de. Substituting for A''^'^ using equation (68) and transforming 
the variable of integration from z to t using equation (24) yields the alternative form 

/■oo 

$f (to, eo, e) = nrlca^ / t^/g(to, t, eo, e) dr , (72) 

Jo 

where is computed using equation (59) and we have also employed equation (26). Despite 
the appearance of the factor in the integrand in equation (72), the contribution to the 
integral from large values of r is actually negligible because the spectrum declines exponen- 
tially in the upstream region due to advection. Furthermore, the escaping spectrum is also 
strongly attenuated in the downstream region (r — > 0) due to the divergence of the electron 
density. Consequently most of the radiation is emitted from the column around r ~ 1 (see 
Fig. 11 from Becker 1998). 

Equations (59) and (72) can be combined to show that the closed-form solution for the 
column-integrated Green's function is given by 



$f (ro, eo, e) = 3iVo kT, e^'^^o ^{eo-e)/{2kn) 

_„ r(i + 2rtr(n + i/2) »"('"°>*^-HH;j""HH7j • '"^ 



where emin, Cmax, K, and /i are defined by equations (40), (47), and (60), and we have 
made the definition 

/■oo 

Xn = / T^gnir) dr . (74) 



The integral X„ can be evaluated analytically by substituting for Qn using equation (44) and 
employing equation (7.414.7) from Gradshteyn & Ryzhik (1980) and equation (15.3.3) from 
Abramowitz & Stegun (1970). After some algebra, the result obtained is 

_ 2 r(n + 1/2) (3 - w)^-\3 -w- Anw) 
" - n! a3/2(3 + w)-+3/2 ' ^ ' 

In Figure 4 we use equation (73) to plot the dependence of the column-integrated Green's 
function on the radiation energy e using the same parameter values employed in Figure 3, 
along with several values for the source optical depth tq, where the monochromatic seed 
photons are injected. In each case, we set a\\ — 10~^a^, a± — a^, a — 0.1 cr^, A^o — 
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Fig. 4. — Column-integrated Green's function $p (s~^ keV~^) describing the photon number 
spectrum escaping from an accretion column as a function of the radiation energy e (keV) 
and the source optical depth tq, evaluated using eq. (73). The curves were computed by 
setting cr|| = 10^^cr.j,, a± = a^, a = 0.1 a.^, Nq = 1, ^ = 2/^/3, tq = 10*^ cm, and = 10'^ K. 
The values of tq, eo, and a are indicated for each plot. 

1, ,^ = To = 10^ cm, and Tg = 10^ K. For the cases with eo = 0.1 keV, thermal 

Comptonization produces a Wien hump at e ~ 1 keV, which is the same behavior noted 
in Figure 3. When a = 0.4, the Wien hump is less prominent due to the effect of bulk 
Comptonization. In general, smaller values of the source optical depth tq result in more 
thermal Comptonization because the scattering plasma is more dense. Conversely, for large 
values of tq, bulk Comptonization tends to produce a power-law spectral shape around the 
photon injection energy. In each case, numerical integration of •I'f with respect to the 
radiation energy e confirms that the total number of photons escaping from the column per 
unit time is equal to unity, which is correct since we have set A^^o = 1. 
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5.4. Column-Integrated Spectrum for an Arbitrary Source 

The vertically-integrated photon spectrum emitted through the walls of the accretion 
column due to an arbitrary source Q is given by (cf. eq. [71]) 

$,(e)= / N,{z,e)dz, (76) 

where TV^ is the altitude-dependent particular solution for the escaping photon spectrum 
computed using equation (70). By substituting for TV^ using equation (70), interchanging 
the order of integration, and applying equation (71), we find that the particular solution for 
the column-integrated spectrum can be written as 

POO POO (f)G(^ f f) 

^e(6) - / / ' ^ °; °' ' el Q{zo, eo) deo dzo , (77) 
Jo Jo iVo 

where $^ is evaluated using equation (73), with tq computed in terms of zq using equa- 
tion (28). Equation (77) facilitates the calculation of the total photon spectrum emitted by 
the entire column as a function of the photon energy e, for any source Q. In § 7 we examine 
the nature of the source term for the various emission processes important in X-ray pul- 
sar accretion columns, and we derive associated results for the column-integrated escaping 
photon spectrum, Astrophysical applications of these results and comparisons with the 
X-ray data for specific pulsars are presented in § 8. 



6. MODEL PARAMETERS AND CONSTRAINTS 

Our model includes a number of free parameters, such as the column radius, tq, the 
temperature of the gas in the thermal mound, 7th, the temperature of the electrons in the 
optically thin region above the mound, Tg, the magnetic field strength, B, and the accretion 
rate, M. There are also several additional theory parameters that are related to the above 
mentioned physical parameters, including a, S, ^, and the scattering cross sections a±, a\\, 
and a. In this section we investigate the relationships between these various quantities, 
and we also introduce several new constraints that are used to reduce the number of free 
parameters. Based on these results, we show that for given values of the stellar mass 
and radius R^, the only quantities that need to be varied when fitting the spectral data for 
a particular source are the six parameters Tg, M, Tq, -B, S, and ^. 
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6.1. Dynamical Constraints 

At large distances from the star, where radiation pressure effects become negligible, we 
require that the approximate velocity given by equation (29) equal the local free-fall velocity. 
We use this condition to calculate the altitude, -Zmax, at the upper surface of the radiating 
region within the accretion column by writing 

2GM, V'"^ 



p , CaTmax , (78) 

where the left-hand side represents the free-fall velocity from infinity to the top of the 
radiative zone. The corresponding optical depth, Tmax, is related to Zmax via (see eq. [28]) 

By using equation (79) to substitute for Tmax in equation (78), we are able to derive a 
quadratic equation for z^bdc with solution 



-S-max — 2 



l + Ci)'/'-l , (80) 



where 



(81) 

We will use this result to set the upper limit for the spatial integrations performed in § 7 when 
we calculate the emergent spectra resulting from the Comptonization of bremsstrahlung and 
cyclotron seed photons. 



6.2. Scattering Cross Sections 

Several electron scattering cross sections are incorporated into the model, based on the 
direction of propagation of the photon. The cross sections for photons propagating parallel 
and perpendicular to the magnetic field direction are denoted by and a±, respectively, 
and the angle-averaged cross section is a. Following Wang & Frank (1981), we set the 
perpendicular cross section equal to the Thomson value, so that (see eq. [7]) 

a±^a^ . (82) 

For given values of the parameters M, tq, and ^, we can compute the parallel cross section 
(T|| by using equation (26) to write 

wrpmpc y 1 
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Once values are specified for the electron temperature and the parameters a and 5, we 
can compute the angle-averaged cross section a by using equation (38), which yields 

"^Sl^T' ^^^^ 

where a\\ is evaluated using equation (83). We shall use equations (82), (83), and (84) to 
set the values of the three scattering cross sections appearing in our model. In general, we 
expect that a\\ -C a -C a± (e.g.. Canuto et al. 1971), and this result will be verified when 
specific numerical models are developed. We will also confirm that the values obtained for 
(T|| are reasonably consistent with the dependence on the mean photon energy expressed by 
equation (6). 



6.3. Thermal Mound Properties 

The blackbody surface of the thermal mound represents the effective photosphere for 
photon creation and destruction in the accretion column, which occurs mainly via free-free 
emission and absorption. The altitude at the top of the mound can therefore be estimated 
by setting the free-free optical thickness across the column equal to unity, 

T^ = ro4 = l, (85) 

where denotes the Rosseland mean of the free-free absorption coefficient. In principle, 
we should use an expression for q;| that accounts for the effect of the strong magnetic field 
(e.g., Lauer et al. 1983). However, the modifications introduced by the presence of the 

field arc relatively minor for photons with energies e ^ kT^, which are the most strongly 
absorbed. We can therefore estimate the value of at the mound surface in the case of 
pure, fully-ionized hydrogen by using equation (5.20) from Rybicki & Lightman (1979) to 
write 

al = 6.10 X 10^^ T,l'/' pI (86) 

in cgs units, where Tth and pth represent the temperature and density at the top of the 
mound, respectively, and wc have set the Gaunt factor equal to unity. 

By combining equations (85) and (86), we find that the density at the thermal mound 
surface is given by 

Ah = 4.05 X 10-^2 Tli' r-'/' . (87) 

The inflow speed at the mound surface, = — '^(^th), is related to pth via the continuity 
equation (19), which yields 

^;th = 7.86 x 10^° Mr^^^^ T~^^^ . (88) 
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We can now compute the optical depth at the top of the thermal mound by combining 
equations (29), (32), and (88), which yields 



rth = 2.64 X 10^« „ ^3/2;7/4. ■ (89) 



MR, 

The corresponding altitude, Zth, is obtained by combining equations (28), (32), and (89), 
which yields in cgs units 

^th = 5.44 X 10^^ . (90) 

We will use Zth as the lower limit for the spatial integrations performed in § 7 when we 
calculate the emergent spectra resulting from the Comptonization of bremsstrahlung and 
cyclotron seed photons. 

The temperature of the gas in the thermal mound, Tth, is computed separately from the 
temperature of the electrons, Tg, in the optically thin region above the mound, although the 
two temperatures are expected to be comparable. We calculate T^h by setting the downward 
flux of kinetic energy at the mound surface equal to the upward flux of blackbody radiation, 
i.e., 

'^s.Ti^^Jvl, (91) 

where 

M 

J = Pth vth = — 2 (92) 

denotes the (constant) mass accretion flux and is the Stephan-Boltzmann constant. 
The right-hand side of equation (91) represents the flux of kinetic energy that must be 
thermalized and radiated by the mound before the gas can merge onto the star. By combining 
equations (19), (88), and (91), we find that the mound temperature is given in cgs units by 

Tth = 2.32 X 10^ M^/s To . (93) 

Once values for the fundamental free parameters M*, i?*, M, ro, and ^ have been selected, 
the mound temperature Tth is evaluated using equation (93), after which the velocity Vt}^, 
optical depth Tth, and altitude Zth at the mound surface are calculated using equations (88), 
(89), and (90), respectively. 



6.4. Bulk vs. Thermal Comptonization 



The dimensionless parameter 5 defined by equation (38) plays a central role in deter- 
mining the overall shape of the reprocessed radiation spectrum resulting from the bulk and 
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thermal Comptonization of the seed photons. We demonstrate here that 5 is essentially the 
ratio of the "|/-parameters" for these two processes. In the one-dimensional flow configura- 
tion under consideration here, the mean energization rate for a photon with energy e due to 
the bulk compression of the decelerating electrons is given by (e.g., Becker 1992) 



de 
di 



_ e dv _ acnecr\\ e 
bulk 3 dz 3 ' 



where the final result is obtained by utilizing equation (24), along with the approximate 
velocity profile given by equation (27). The associated y-parameter, describing the average 
fractional energy change experienced by a photon before it escapes through the column walls, 
is defined by 

^esc / ds 



bulk 3 

where tesc is the mean escape timescale, expressed as a function of z by equation (65). 



(95) 



Next we recall that the mean stochastic energization rate for the thermal Comptoniza- 
tion of photons is given by (Rybicki & Lightman 1979) 



de 
dt 



^Ueac e ^ , (96) 

therm "^eC 



with associated |/-parameter defined by 



^esc / da 



ytheriii — \ 7/ 

dt 



= rie (7 C tesc ^ ■ (97) 

therm mgC^ 



We note that Ue oc v"^ oc z"^^"^ and tesc z"^/^ according to equations (29) and (65), 
respectively. Hence, ybuik and ytherm each vary in proportion to z~^, and therefore both 
quantities diverge at the base of the column. It is therefore difficult to assign global values 
for the two parameters that characterize the entire accretion column. 

We can gain insight into the relative importance of bulk and thermal Comptonization by 
relating the ratio of the respective y-parameters to the value of 5 by combining equations (38), 
(95), and (96), which yields 

^ = ^ . (98) 

4 1/thcrm 

According to this relation, bulk Comptonization dominates the photon energization if 5 ^ 1, 
and thermal Comptonization dominates if 5 <C 1. In the intermediate regime, with 5 ~ 1, 
both processes contribute about equally to the energy transfer between the electrons and 
the photons. However, it is important to realize that even for large values of 5, the thermal 
process may still have a profound effect on the overall shape of the spectrum by transferring 



energy from high to low frequency photons, which produces a quasi-exponential cutoff at 
high energies and a concomitant flattening of the spectrum at lower energies. 



By combining equations (94) and (96) with equation (7.36) from Rybicki & Lightman 
(1979), we flnd that the total (mean) photon energization rate for the thermal+bulk Comp- 
tonization model considered here is given by 

' de 



dt/ total \^eC^ rUeC^ 3(7/ 

where the terms on the right-hand side describe the effects of stochastic energization, electron 
recoil losses, and bulk Comptonization, respectively. The total energization rate vanishes at 
the critical energy 

eerit = 4A;re + ^mec2 . (100) 

O (J 

In the low-temperature limit, we flnd that bulk Comptonization balances recoil losses at the 
energy 

ebuik=^mec% Te^O, (101) 

which is the maximum photon energy achieved in the thermal-|-bulk model as Tg — > 0. Note 
that in the pure bulk Comptonization model of Becker & Wolff (2005a,b), recoil losses are 
not included, and therefore the power-law shape extends to inflnite photon energy. This 
makes it necessary to restrict the index of the photon spectrum to ax > 2 in order to avoid 
an infinite photon energy density in the Becker & Wolff model. However, no such restriction 
exists in the thermal+bulk model considered here because recoil losses always attenuate the 
spectrum at high photon energies. 



6.5. Photon Advection, Diffusion, and Escape 

The seed photons injected into the accretion column experience a variety of physical 
effects as they propagate within the plasma, eventually escaping through the column walls. 
In this section we will show that in a radiation-dominated accretion column, the timescale for 
the escape of radiation through the column walls is comparable to the dynamical timescale 
for the accretion of the gas onto the star. In particular, we will demonstrate that the 
dimensionless parameter ^ defined in equation (26) is proportional to the ratio of these two 
timescales. 

We define the dynamical timescale, tshock, for the accretion of material from the sonic 
point down to the surface of the star by writing 
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where Zsp is the distance between the sonic point and the stellar surface given by equa- 
tion (31). We can form the ratio of the dynamical timescale tghock to the mean escape 
timescale tesc for photons to diffuse through the walls of the accretion column by employing 
equations (18) and (102), which yields 



^shock ^ C Zgp 

tesc M (7 ± 



(103) 



We can substitute for M and Zgp using equations (26) and (31), respectively, to obtain the 
alternative result 

W^ln(7y3) ^^Q_24e. (104) 

^esc 2 V 3 

Becker (1998) found that ^ = 2/v^ in a radiation-dominated pulsar accretion flow, and 
therefore we conclude that 

ishock ~ 0.28 iesc ■ (105) 

This result is expected, since the timescale for the photons to escape from the column must 
be comparable to the dynamical timescale because the radiation must be able to remove the 
kinetic energy from the flow. This is a general property of radiative accretion flows (see, 
e.g., Imamura et al. (1987). 

We can gain additional insight into the processes affecting the transport of photons 
through the accretion column by analyzing the relative importance of the upward diffusion 
of photons along the column axis, compared with the downward advcction of photons caused 
by the "trapping" of radiation in the rapidly falling gas. The parameter ^ defined in equa- 
tion (26) also provides a direct measurement of the inffuence of these two competing effects. 
The mean diffusion velocity for photons propagating parallel to the axis of the accretion 
column, w\\, can be approximated by writing 

W|| = ^ . (106) 

Photon "trapping" occurs when the downward advective ffux dominates over the upward 
diffusive flux, so that the photons are essentially conflned in the lower region of the flow 
(Becker & Begelman 1986). In the model considered here, we flnd that the radiation is 
trapped if \v\ > w\\, which can be rewritten using equations (27) and (106) as the condition 
r > Ttrap, where 



In the lower region of the column, where r < Ttrap, diffusion is able to transport radiation 
effectively in the vertical direction. We can calculate the associated "trapping altitude" for 
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photons diffusing vertically in the column by combining equations (28) and (107), which 
yields 

..„^ = !Ill£!!!. = l5i('^V'\ (108) 



2M(7|| 

where the final result is obtained by substituting for M using equation (26). It is interesting 
to compare the value of 2;trap with the distance from the stellar surface to the sonic point, 
Zs^, given by equation (31). We find that the ratio of these two quantities is given by 

^-r^-2.oe. (109) 

ln(7/3) 

Since ^ ~ 1 in a radiation-dominated pulsar accretion flow, equation (109) indicates that z^^ 
and 2;trap are comparable, and therefore the photons are not able to penetrate very far into 
the region above the sonic point. Hence most of the photons escape through the walls of the 
column below the sonic point, with a mean residence time equal to the dynamical time for 
the gas to settle onto the stellar surface (see eq. [105]). 



6.6. Radius of the Accretion Column 

In our model, the radius Tq of the accretion column (which is also the radius of the "hot 
spot" on the stellar surface) is treated as a free parameter. However, it is possible to constrain 
the value of this radius based on dynamical considerations related to the magnetically- 
channeled accretion of the gas onto the neutron star. For example, by taking into account 
the effects of both the stellar rotation and the strong magnetic held. Lamb et al. (1973) find 
that the radius of the hot spot satisfies the constraint (see also Harding et al. 1984) 

/D \i/2 

ro < i?J — , (110) 



where is the stellar radius and ta is the Alfven radius, given in cgs units by 

-■A-Z-exlO (jjlJcj (^j (^j ( io^Terg.B-0 ■ 

By combining equations (9), (110), and (111), we find that column radius tq is bounded by 
the upper limit 



1012 Gy VlOkm/ \MqJ Il0i7gs-i 

which represents an interesting constraint for our model. We will compare our results for tq 
with this upper bound when specific numerical models are developed in § 8. 
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6.7. Compcirison of Exact and Approximate Velocity Profiles 

Following Lyubarskii & Sunyaev (1982), we have adopted the approximate velocity 
profile given by equation (29), which faciliates the separation of the transport equation. The 
approximate velocity profile was compared with the exact solution (eq. [30]) in Figure 2, 
from which we conclude that the two functions agree fairly well in the radiative portion of 
the accretion column. We can further explore the validity of equation (29) by comparing the 
results obtained for the emergent spectra using the exact and approximate velocity profiles. 
This can be accomplished by contrasting the spectra computed using the thermal+bulk 
Comptonization model developed here (based on the approximate velocity distribution) with 
tho corresponding results obtained using the "pure" bulk Comptonization model discussed 
by Becker & Wolff (2005b), which is based on the exact velocity profile. 

In Figure 5 we display the results obtained for the column- integrated Green's function, 
calculated using either the thermal+bulk Comptonization model (eq. [73]) or the pure 
bulk Comptonization model (eq. [81] from Becker & Wolff 2005b). The parameters values 
employed are the same ones adopted in Figure 3, with Cq = 0.01 keV, a = 0.4, ay = 10^^ o"^, 
a± = a^,a = 0.1 a^, Nq = 1, ^ = 2/\/3, tq = 0.1, and Tq = W cm. The electron temperature 
Tg is set equal to 10^ K, 10^-^ K, 10^-^ K, or 0, and the corresponding values for S are 0.79, 
3.95, 39.54, and oo, respectively (see eq. [38]). Note that a Wien hump forms at energy 
e ~ ^crit (eq. [100]) for large values of Tg due to photon energy redistribution, and as Tg — > 
the spectrum approaches the power-law shape characteristic of bulk Comptonization. The 
high-energy cutoffs displayed by the spectra in Figure 5 are the result of electron recoil only; 
in actual X-ray pulsar spectra, the quasi-exponential cutoffs are produced by a combination 
of electron recoil and cyclotron absorption, as discussed in § 8. 

The Becker & Wolff (2005b) spectrum is based on the exact velocity solution (eq. [30]), 
and therefore it provides an interesting test for the model considered here, which is based on 
the approximate velocity profile (eq. [29]). In the limit Tg — > 0, bulk Comptonization domi- 
nates the spectral formation process, and the column-integrated Green's function computed 
using equation (73) agrees perfectly with the Becker & Wolff result up to photon energy 
e = Cbuik, defined as the energy at which bulk Comptonization is balanced by losses due to 
electron recoil (see eq. [101]). Since recoil losses were not included in the Becker & Wolff 
(2005b) model, the resulting spectrum extends to infinite photon energy. Similar agreement 
between the two models in the low-temerature limit has been confirmed using several other 
sets of parameters. These comparisons confirm that the approximate velocity profile uti- 
lized here provides a good representation of the effect of bulk Comptonization in the pulsar 
accretion column. 
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Fig. 5. — Column-integrated Green's function $p (s~^ keV~^) describing the photon number 
spectrum escaping from an accretion column as a function of the radiation energy e (keV) 
and the electron temperature T^, evaluated using eq. (73). The results were computed by 
setting eo = 0.01 keV, a = 0.4, ct|| = 10^^ a^, a±_ = a^, a = 0.1 a^, iVo = 1, ^ = 2/-\/3, 
To = 0.1, and vq = 10"^ cm. The corresponding value of Tg is indicated for each curve, and 
the heavy dashed line denotes the spectrum obtained using the exact velocity profile; see the 
discussion in the text. 

7. PHOTON SOURCES AND ASSOCIATED SPECTRA 

The high-energy radiation spectrum emerging from an X-ray pulsar accretion column is 
produced via the bulk and thermal Comptonization of seed photons injected into the plasma 
by a variety of emission mechanisms. Based on the linearity of the transport equation (15), we 
can compute the radiation spectrum emerging from the accretion column for any source term 
Q by employing the integral convolutions given by equations (70) and (77) for the altitude- 
dependent and column-integrated spectra iVe(z,e) and $e(e), respectively. The calculation 
of the two spectra is straightforward due to the availability of the analytical solutions for 
the altitude-dependent Green's function A*"^^ (eq. [68]) and for the column-integrated Green's 
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function (eq. [73]). In the context of accretion-powered X-ray pulsars, the primary sources 
of seed photons are bremsstrahlung, cyclotron, and blackbody radiation (Arons et al. 1987). 
The first two mechanisms inject photons throughout the column, although the emission tends 
to be concentrated towards the base of the column because of the strong density dependence 
of these processes. Conversely, blackbody radiation is injected primarily at the surface of the 
thermal mound (i.e., the "photosphere"), where the gas achieves thermodynamic equilibrium. 
The energy dependences of the three emission mechanisms are also quite different, with the 
bremsstrahlung and blackbody processes producing broadband continuum radiation, and the 
cyclotron producing nearly monochromatic radiation. 

The nature of the source term Q appearing in the transport equation (15) is examined 
in detail below for each of the emission mechanisms of interest here. In general, Q can be 
related to the photon emissivity, rig, by writing 



where de expresses the number of photons created per unit time per unit volume in the 
energy range between e and e + de. We will use equation (113) to determine Q by specifying 
the emissivity fie for bremsstrahlung, cyclotron, and blackbody emission. In our approach, we 
treat the cyclotron and bremsstrahlung sources as separate emission mechanisms, although 
they can also be viewed as components of a single, generalized magneto-bremsstrahlung 
process, as discussed by Riffert et al. (1999). 



In an accretion-powered X-ray pulsar, cyclotron photons are produced mainly as a result 
of the collisional excitation of electrons to the n — 1 Landau level, followed by radiative 
decay to the ground state (n = 0). The electrons are excited via collisions with protons, and 
therefore the production of the cyclotron radiation is a two-body process. Estimates suggest 
that in the luminous sources, collisional dccxcitation is completely negligible compared with 
radiative deexcitation, and therefore the production of the cyclotron radiation acts as a 
coohng mechanism for the gas (e.g., Nagel 1980). The subsequent Comptonization of these 
photons results in further coohng. The cyclotron hue will be broadened by thermal effects, 
as well as the possible spatial variation of the magnetic field, but these effects are likely to be 
masked by the strong diffusion in energy space caused by thermal Comptonization. We can 
therefore approximate the distribution of the cyclotron seed photons using a monochromatic 
source term. Based on equations (7) and (11) from Arons et al. (1987), we find that in the 
case of pure, fully-ionized hydrogen, the rate of production of cyclotron photons per unit 




(113) 



7.1. Cyclotron Radiation 
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volume per unit energy is given in cgs units by 



nf" = 2.10 X 10^^ p2 ^-3/2 ^ 




) 



e 



(114) 



where the cyclotron energy Cc is defined by equation (1), and 



H 



( 




0.41 , ec/kTe > 7.5 , 



0.15 y/eJkTe , eJkTeK 7.5 . 



(115) 



The cyclotron fine will be broadened by thermal effects, as well as the possible spatial 
variation of the magnetic field, but these effects are hkely to be masked by the strong 
diffusion in energy space caused by thermal Comptonization. We can therefore approximate 
the cyclotron source term, Q^^^, using the monochromatic expression 



obtained by combining equations (113) and (114). 

The cyclotron source term described by equation (116) is localized in energy and dis- 
tributed in space. Following Lyubarskii & Sunyaev (1982), we assume that the electrons 
in the optically thin plasma above the thermal mound are isothermal, and we also assume 
that the magnetic field B has a cylindrical geometry with a constant value throughout the 
column. The neglect of the dipole variation of the field geometry and strength is expected 
to have a small effect on the emergent spectrum because most of the radiation is created, re- 
processed, and emitted through the column walls within one or two scattering optical depths 
above the stellar surface. Compared with the cylindrical geometry employed in our model, 
the inclusion of the actual dipole geometry of the magnetic field will tend to increase the 
plasma density at the bottom of the accretion column, due to the funnel-like shape of the 
column walls. This is likely to further increase the concentration of the radiative processes 
near the bottom of the column. However, a small fraction of the photons may nonetheless 
diffuse to a significant altitude (perhaps several kilometers) above the star before escaping, 
because the scattering cross section for photons propagating parallel to the magnetic field 
is generally much less than Thomson (see § 6.2). This issue requires further exploration 
using a numerical model that includes the exact dipole dependence of the magnetic field. 

For a constant value of B, the spatial variation of the cyclotron source term reduces to 
Q^^^ oc p^, which reflects the two-body nature of the coUisional excitation process. We can 
obtain the particular solution for the column-integrated spectrum of the escaping radiation 
resulting from the reprocessing of the cyclotron photons by using equation (116) to substitute 
for Q in equation (77). The energy integration is trivial, and we find after substituting for 



Q'^'iz, e) = 1.92 X 10^2 ^2 ^2 ^-7/2 ^ 




■■"^^ S{e - 6e) , 



(116) 



-41 - 



Bi2 using equation (1) and integrating over zq that the particular solution for the column- 
integrated escaping photon spectrum is given in cgs units by 



^+3/2 {e,+e)/(2fcTe: 



n=0 



r(/i- K + 1/2) n! 
r(l + 2/i)r(n+l/2) 



X X„ Ar, iV4,^[min(x, Xc)] W^^^[m&x{x: Xc)] , 



(117) 



where X„ is computed using equation (75) and we have made the definitions x = e/ZcTg and 
Xc = ^c/kTf.. The quantity A„ represents the spatial integration 



An= I e 

Jrth 



3ar|/2 



5'n(To) To dro , 



(118) 



where Tmax and Tth denote the optical depths at the top of the accretion column and the 
surface of the thermal mound, respectively (see eqs. [79] and [89]). By substituting for gn 
using equation (44), we find that An can be evaluated analj^ically, yielding 



_ -A r(n + 1/2) (2^/)™ (3 - w)- 

An- 



r(m+ l/2)2m!(n-m)! 



m. 



a{w - 3)Tti 



a{w — 3)t, 



2 

max 



(119) 



where r(a, z) denotes the incomplete gamma function. 



7.2. Blackbody Radiation 

Photons are produced at the surface of the thermal mound with a blackbody distribu- 
tion, and therefore the surface flux is equal to tt times the intensity (Rybicki & Lightman 
1979). Following Becker & Wolff (2005b), we deflne the function S{e) so that e^S{e)de 
represents the number of photons emitted from the surface per second in the energy range 
between e and e+de. We can relate S{e) to the Planck distribution by expressing the amount 
of energy emitted per second from the surface of the thermal mound (with area nvo) in the 
energy range between e and e + de as 

e^S{e)de^TTrl-nB,{e)de , (120) 

where 

denotes the blackbody intensity and Tth is the temperature of the gas in the thermal mound. 
Note that the units for 5^ are ergs s~^ ster~^ cm~^ erg"-^. 
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The function S is related to the source term Q appearing in the transport equation (15) 



via 



Q''''{z,e)=S{e)S{z-z,i,) , (122) 
which can be combined with equations (120) and (121) to conclude that 



This result for the blackbody source term is localized in physical space, but it possesses a 
distributed (continuum) energy dependence, which is the exact opposite of the cyclotron 
source given by equation (116). 

We can obtain the particular solution for the column-integrated spectrum of the escaping 
radiation resulting from the reprocessing of the blackbody photons by using equation (123) 
to substitute for Q in equation (77). In this case the spatial integration is trivial, and we 
find after some algebra that 

chbb/.N _ 2 7rVg ^f(^th,eo,e) egdep 



c'^h^ Jo No e^o/'^^th _ 1 ' 

where $^ is computed using equation (73). 



7.3. Bremsstrahlung Radiation 

The bremsstrahlung emission spectrum produced by electrons streaming along a strong 
magnetic field is modified from the standard thermal result due to the appearance of the reso- 
nance at the cyclotron energy Cc (Lauer et al. 1983). However, it can be shown that the angle- 
averaged emission spectrum obtained by integrating over the onc-dimcnsional Maxwellian 
electron distribution in the accretion column agrees fairly well with the field-free thermal 
expression for photon energies e < Cc- We can therefore describe the bremsstrahlung pho- 
ton source using the standard thermal formula in combination with a separate (resonant) 
cyclotron term described by equation (114). Based on equation (5.14b) from Rybicki & 
Lightman (1979), we find that the thermal component of the brcmsstrahhmg (free-free) 
photon production rate per unit volume per unit energy is given in cgs units by 

nf = 3.68 X 10^6 p2 T-V2 e'^ e'^/'^^^ , (125) 

for a plasma composed of pure, fully-ionized hydrogen. Equation (125) can be combined 
with equation (113) to define the free-free source term 

Q''{z,e) = 1.16 X 10=^VoV'7;-V2e-3e-^A^« , z^^ < z < w , (126) 
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where Zmax and Zth denote the altitudes at the top of the accretion column and the surface 
of the thermal mound, respectively (see § 6.1). Equation (126) describes the spectrum of 
the injected, optically-thin bremsstrahlung emission for photons with energy eabs < e < ec, 
where eabs represents the thermal self-absorption cutoff, below which the spectrum becomes 
optically thick. The calculation of tabs is complicated by the fact that the soft photons 
created via bremsstrahlung emission may be upscattered due to collisions with fast electrons 
before they are reabsorbed. The energy amphfication associated with electron scattering 
tends to enhance the survival probability of the photons because the absorption length is a 
rapidly increasing function of the photon energy. The effective value of eabs can therefore be 
estimated by setting 4c — l/ctfj where 4c — (^T^^e)"^ is the electron scattering mean free 
path and is the free-free absorption coefficient given by equation (5.18b) from Rybicki & 
Lightman (1979). This yields 

^ = 6.08 X 10^2^-7/4^1/2 _ (^27) 

The utilization of the Thomson cross section is justified because most of the bremsstrahlung 
photons are emitted roughly perpendicular to the magnetic field due to the dipole nature of 
the radiation pattern (e.g., Lauer et al. 1983). In our applications, the value adopted for 
eabs is obtained by averaging equation (127) over the vertical structure of the column. We 
typically find that 0.01 < eabs/^^e ~ 0.1. 

The bremsstrahlung source term given by equation (126) possesses distributed spatial 
and energy dependences, and therefore it is more challenging to implement than the cy- 
clotron and blackbody sources given by equations (116) and (123), respectively. By using 
equation (126) to substitute for Q in equation (77), we find after several steps that the par- 
ticular solution for the column-integrated spectrum of the escaping radiation resulting from 
the reprocessing of the bremsstrahlung photons can be written in cgs units as 

ff, , 2.80 X IQ-^^Mi^V^ g«-2 g-6/(2fcTe) 
$ le) = ^ 

^ r(l + 2/i)r(n + l/2) " " ^ > 



X 



n=0 



where X„ and An are computed using equations (75) and (119), respectively. The quantity 
Bn represents the energy integral 

/■oo 

Bn= Xo M,,^[min(x, Xo)] W^.,^[max(x, Xo)] cixo , (129) 

where x = ^/kT^i Xo = ^o/kT^, and Xabs = ^abs/kTg. This integral can be evaluated analyti- 
cally, but the result is very complex and therefore we do not state it here. 
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8. ASTROPHYSICAL APPLICATIONS 

By combining results from the previous sections in the paper, we can now compute 
the spectrum emitted from an X-ray pulsar accretion column due to the bulk and thermal 
Comptonization of bremsstrahlung, cyclotron, and blackbody seed photons. The theoretical 
prediction for the phase-averaged photon count rate spectrum measured at Earth is computed 
using 

where 

$f (e) = + $bb(^) + $ff(^)] (131) 

denotes the total spectrum of the Comptonized radiation escaping from the column, D is 
the distance to the source, and the function represents a Gaussian cyclotron absorption 
feature given by the form (e.g., Heindl & Chakrabarty 1999; Orlandini et al. 1998; Soong et 
al. 1990) 

A^{e) = 1 - e-(e-ec)V(2<^?) . (132) 

The quantities on the right-hand side in equation (131) express the contributions to the 
escaping spectrum due to the reprocessing of cyclotron radiation (^^y'^; eq. [117]), blackbody 
radiation ($^^; eq. [124]), and bremsstrahlung emission ($f ; eq. [128]), respectively. In this 
section we will use equation (130) to compute observed spectra based on parameters corre- 
sponding to three specific X-ray pulsars. The theoretical spectra will be compared with the 
observational data for each source in order to gain insight into the effects of bulk and ther- 
mal Comptonization in the accretion shock on the formation of the X-ray continuum. The 
spectral results presented here are example calculations, rather than detailed quantitative 
fits, which will be reported in subsequent papers. 

The sources treated here are Her X-1, LMC X-4, and Cen X-3, which represent several 
of the brightest known X-ray pulsars. In each case the values adopted for the stellar mass 
and radius are M* = 1.4 M0 and — 10 km, respectively. The accretion dynamics in these 
high-luminosity sources is expected to be dominated by radiation pressure, and therefore 
we expect that ^ ~ 2/a/3 (see Becker 1998). The scattering cross section for photons 
propagating perpendicular to the magnetic field in our calculations is given by a± = o"^ 
(see eq. [82]), the accretion rate M is constrained by the observed X-ray luminosity, and 
the source distance D is determined using published estimates. The remaining fundamental 
input parameters that must be set in order to calculate the theoretical spectrum are S, B, 
To, and Tg, which are hsted in Table 1 for each of the models considered here. Once these 
quantities are specified, the additional parameters a, Tmax, cr\\, ^, Tth, ^th, ^trap) ^-nd eabs can 
be computed using equations (33), (79), (83), (84), (89), (93), (107), and (127), respectively. 
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with the results reported in Table 2. Note that in each case, we find that Tmax ~ Ttrap, which 
indicates that the observed emission is produced in the "trapped" region of the accretion 
column, as expected (see § 6.5). Table 2 also includes results for the inflow speed at the 
thermal mound surface, Vth (eq. [88]), and the accretion mass flux, J (eq. [92]). In general, 
we find in our applications that fth ~ 0.03 c, which represents an enormous deceleration from 
the free-fall velocity vs ~ 0.6 c at the top of the column. Since there is not much variation 
in the value of I'th, we expect that the mound temperature Tth will increase with increasing 
J (see eq. [91]). 

Once all of the parameters are specified, the corresponding count-rate spectrum is eval- 
uated using equation (130). The spectral results presented here were computed using the 
first five terms in the expansions (eqs. [117], [124], and [128]), which generally provide at 
least three decimal digits of accuracy. As a check on the numerical results for the spectral 
components, we confirm that the number of photons escaping from the column per unit time 
is exactly equal to the number injected, as required in the steady-state scenario considered 
here. In addition to the Comptonized spectrum evaluated using equation (130), the results 
for the total spectra displayed in Figures 6, 7, and 8 also include an Fe Ka emission line 
modeled using the Gaussian function 

F^{e) ^ e-{.-.K)V(2.|) , (133) 

(7kv 27r 



where is the associated total photon fiux measured at Earth. Note that in the approach 
taken here, this spectral feature is treated in an ad hoc manner by simply adding it to the 
final spectrum rather than by subjecting the iron line photons to subsequent Comptonization. 
Unfortunately, a more sophisticated approach is beyond the scope of the present paper due 
to the complexity of the spectral formation process for the iron emission line photons. The 
associated auxiliary parameters ex, (ix, and ctk are listed in Table 3, along with Cc, dc, and 
(7c (see eq. [132]). 



8.1. Her X-1 

In Figure 6 we display the theoretical count-rate spectrum computed using equation (130) 
along with the deconvolved (incident), phase- averaged BeppoSAX data reported by Dal Fl- 
ume et al. (1998) for Her X-1, which has the instrumental response removed. The values 
used for the fundamental theory parameters correspond to model 1, with ^ = 1.45, aj^ = a^, 
S = 1.80, B = 3.80 X 10^2 G, ro = 44m, = 6.25 x 10^ K, M = 1.11 x 10^^gs-^ and 
D — 5kpc, as indicated in Table 1. The associated values for the computed parameters a, 
Tmax, cr\\, ^, Tth, ^ud Tth ^rc reported in Table 2, and the additional auxiliary parameters are 
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listed in Table 3. Results are presented for the total spectrum, as well as for the individual 
contributions to the observed flux due to the Comptonization of cyclotron, blackbody, and 
bremsstrahlung seed photons. The lack of a strong Wien peak in the spectrum, along with 
the moderate value of 6, indicate that thermal Comptonization is unsaturated in this source. 
It is interesting to note that the value wc obtain for tq is about an order of magnitude 
below the upper limit for the column radius computed using equation (112), which yields 
ro < 462 m. 

Although the results presented here arc not fits to the data, it is interesting to note 
that the general shape of the pulsar spectrum predicted by the theory agrees very well with 
the observations for Her X-1, including both the quasi-exponential cutoff at high energies 
and the power law shape at lower energies. In the case of Her X-1, we find that reprocessed 
(Comptonized) blackbody emission from the thermal mound makes a relatively small con- 
tribution to the spectrum due to the small radius of the accretion column. We also find that 
reprocessed cyclotron emission is swamped by reprocessed bremsstrahlung due to the high 
temperature of the infalling plasma. Hence we conclude that the X-ray spectrum emitted 
by Her X-1 is completely dominated by Comptonized bremsstrahlung emission. 



8.2. LMC X-4 

In Figure 7 we compare the theoretical count-rate spectrum computed using equa- 
tion (130) with the deconvolved, phase-averaged BeppoSAX data reported by La Barbera 
et al. (2001) for LMC X-4. The values used for the fundamental theory parameters cor- 
respond to model 2, with ^ = 1.15, = a^, 5 = 1.30, B = 3.28 x lO^^G, tq = 680m, 
= 5.90 X 10^ K, M = 2.00 x 10^^gs-\ and D = 55kpc (see Table 1). The associated 
values for the computed and auxiliary parameters are listed in Tables 2 and 3, respectively. 
Note that the results for the accretion rate M and the column radius ro are much larger 
than those for Her X-1, although the electron temperatures are nearly the same. In this 
case, reprocessed cyclotron emission makes a slightly stronger contribution to the high- 
energy spectrum than in model 1, although the X-ray spectrum is clearly dominated by 
Comptonized bremsstrahlung emission. In model 2, the effect of thermal Comptonization 
(relative to bulk) is somewhat stronger than in model 1, as indicated by the smaller value of 
S (see eq. [97]). This effect tends to flatten the LMC X-4 spectrum in the 2-22 keV energy 
range as compared with the Her X-1 spectrum, although the thermal Comptonization is 
clearly unsaturated, as demonstrated by the absence of a Wien peak in the spectrum. We 
also note that the reprocessed blackbody emission makes a much smaller contribution to the 
observed spectrum in model 2 than in model 1 because the thermal mound temperature 7th 
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Fig. 6. — Theoretical column- integrated count rate spectrum -^^(e) evaluated using eq. (130) 
based on the model 1 parameters listed in Table 1, compared with the deconvolved (incident) 
X-ray spectrum for Her X-1 reported by Dal Fiume et al. (1998; circles and crosses). The 
plots include the total spectrum as well as the individual components due to Comptonized 
bremsstrahlung radiation, Comptonized cyclotron emission, Comptonized blackbody radia- 
tion, and the iron emission line feature, as indicated. Note that the total spectrum is strongly 
dominated by the bremsstrahlung component. 

has dropped substantially, due to the marked decrease in the mass flux, J, as indicated in 
Table 2. In contrast to the Her X-1 case, for LMC X-4 we find that the value obtained for r^ 
is comparable to the upper limit for the column radius computed using equation (112), which 
gives ro < 727 m. This suggests a possible connection between the column radius and the 
luminosity in accretion-powered X-ray pulsars that can be further explored using parameter 
studies based on a large number of sources. 
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Fig. 7. — Same as Fig. 6, except the data correspond to LMC X-4 and the theoretical 
spectra were computed based on the model 2 parameters reported in Table 1. The data were 
reported by La Barbera et al. (2001; circles and crosses). 

8.3. Cen X-3 

In Figure 8 we plot the deconvolved, phase-averaged BeppoSAX data reported by Bur- 
deri et al. (2000) for Cen X-3 along with the theoretical count-rate spectrum computed 
using equation (130). The values adopted here for the fundamental theory parameters cor- 
respond to model 3, with ^ = 1.25, = a^, 6 = 3.71, B = 2.63 x lO^^G, tq = 730m, 
Te = 3.40 X 10^ K, M = 1.51 x 10^^gs-\ and D = 8kpc (see Table 1), and the additional 
parameters are listed in Tables 2 and 3. In this case only, interstellar absorption has been 
included based on a hydrogen column density of A^h = 2.0 x 10^^ cm~^ in order to reproduce 
the observed low-energy turnover (cf. Burderi et al. 2000). We note that the values obtained 
for the accretion rate M and the column radius ro are once again much larger than those for 
Her X-1. It is apparent that the contributions due to Comptonized cyclotron and blackbody 
emission are completely negligible, and the observed spectrum is dominated by Comptonized 
bremsstrahlung photons. 
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Fig. 8. — Same as Fig. 6, except the data correspond to Cen X-3 and the theoretical spectra 
were computed based on the model 3 parameters reported in Table 1. The data were reported 
by Burderi et al. (2000; circles and crosses). This case also includes interstellar absorption 
with hydrogen column density A^h = 2.0 x 10^^ cm~^. 

The effect of bulk Comptonization is more pronounced in this case than in the spectra 
of Her X-1, as indicated by the larger value obtained for 6. We also note that the values 
obtained for the accretion flux J and the thermal mound temperature Tth in model 3 are 
lower than those associated with models 1 and 2. However, despite the decrease in the mound 
temperature, the Comptonized blackbody emission has increased in strength to become 
roughly comparable to the reprocessed cyclotron emission, which reflects the decrease in 
the magnetic field strength B. As in the case of LMC X-4, we find that the value for tq 
associated with Cen X-3 is comparable to the upper limit for the column radius computed 
using equation (112), which yields tq ^ 744 m. This supports the possibility of a connection 
between the column radius and the X-ray luminosity, as suggested in § 8.2. 
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9. CONCLUSIONS 

The new theoretical model developed here represents the first complete calculation of the 
X-ray spectrum associated with the physical accretion scenario first suggested by Davidson 
(1973), including bulk and thermal Comptonization. For given values of the stellar mass 

and the stellar radius i?*, our model has six fundamental free parameters, namely the 
accretion rate M, column radius ro, the electron temperature Tg, the magnetic field strength 
B, the photon diffusion parameter ^ (eq. [26]), and the Comptonization parameter S (eq. [38]). 
The accretion rate is constrained by the observed luminosity, which effectively removes this 
parameter from the set. A unique solution for the remaining free parameters can be found by 
comparing the model spectrum with the observations for a particular source. The analytical 
nature of the new model makes it amenable to incorporation into standard analysis packages 
such as XSPEC. 

We have shown that the combination of bulk and thermal Comptonization naturally 
leads to the development of emergent spectra in accretion-powered X-ray pulsars that are in 
good agreement with the observational data for three bright sources. The main result of the 
derivation is the analytical solution for the column-integrated Green's function describing 
the escaping radiation spectrum, given by equation (73). Based on this closed-form solution, 
we have developed expressions for the emergent spectra resulting from the reprocessing of 
seed photons injected as a result of cyclotron, blackbody, and bremsstrahlung emission, given 
by equations (117), (124), and (128), respectively (we note that additional seed photons may 
also be supplied via external illumination of the column, which is not considered here). In 
each case, we confirm that the number of photons escaping from the column per unit time 
is exactly equal to the number injected, as expected in the steady-state scenario considered 
here. Our main results, summarized in Figures 6, 7, and 8, indicate that the spectra observed 
from the bright sources Her X-1, LMC X-4, and Cen X-3 are dominated by Comptonized 
bremsstrahlung. 

The treatment of the full energy and angle dependence of the scattering cross sections 
for the two polarization modes computed using equations (2) and (4) is not possible using 
an analytical approach. As an alternative, in this paper we have adopted a two-dimensional 
approximation by modeling the scattering of photons propagating either parallel or per- 
pendicular to the magnetic field using the energy- and mode-averaged expressions given by 
equations (6) and (7), respectively. The results we have obtained for the directional compo- 
nents of the cross section listed in Table 2 are not unreasonable, given the uncertainties in 
the detailed physical properties of the accretion column. In particular, we find that the cross 
section for photons propagating parallel to the field direction, (T\\, is substantially reduced 
relative to the Thomson value. This result is consistent with equation (6), since e -C eg for 
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most of the photons in typical X-ray pulsar spectra. Furthermore, the cross section values 
reported in Table 2 confirm that cry <C a -C aj_, as expected (see § 6.2 and Canuto et al. 
1971). Compared to the behavior obtained in our model, we expect that utilization of the 
exact cross sections given by equations (2) and (4) would tend to increase the efficiency of 
thermal Comptonization, relative to the bulk process, for photons with energy e close to the 
cyclotron energy ec- This effect is therefore expected to cause some modification of the pa- 
rameter values required to fit the observed X-ray pulsar spectra. Since the radiative transfer 
problem including the exact cross sections is not tractable analytically, further study of this 
issue awaits the development of a complete numerical model, which we intend to pursue in 
future work. 

The validity of the approximate velocity profile employed here (eq. [29]) was established 
in Figure 5 by demonstrating excellent agreement between the bulk Comptonization spec- 
trum computed using our new model in the limit Te — > and the corresponding spectrum 
obtained using the pure bulk Comptonization model of Becker & Wolff (2005b), which is 
based on the exact velocity profile (eq. [30]). The values obtained for the accretion column 
radius tq (reported in Table 1) comply with the corresponding upper limits computed using 
equation (112), which was derived based on the detailed dynamics of neutron star accretion 
flows (Lamb et al. 1973; Harding et al. 1984). We find that Tq is essentially equal to the 
upper hmit for the two brightest sources considered here (Cen X-3 and LMC X-4), while it 
is about an order of magnitude smaller than this limit for Her X-1. Since the luminosity of 
Her X-1 is roughly an order of magnitude lower then the other two sources, this suggests 
a possible connection between the luminosity and the column size that should be further 
investigated by developing numerical models that simultaneously treat the dynamics and 
the spectral formation in X-ray pulsars. 

The transfer of energy from the gas to the photons is generally dominated by bulk Comp- 
tonization in accretion-powered X-ray pulsars. However, our results confirm that unsaturated 
thermal Comptonization nonetheless plays a significant role in the spectral formation pro- 
cess. Specifically, we find that the thermal process effectively transfers energy from high to 
low frequency radiation, which contributes to the observed quasi-exponential cutoffs at high 
energies and at the same time causes flattening of the spectrum at lower energies. While 
it has been known for some time that pulsars generally have spectra that are well fitted 
using a combination of a power-law plus a quasi-exponential cutoff, until now this form has 
been adopted in a purely ad hoc manner. The new model described here finally provides 
a firm theoretical foundation for this empirical result, based on model parameters that are 
directly tied to the physical properties of the source. It therefore represents a significant 
step in the development of a comprehensive theory for the spectral formation process in 
accretion-powered X-ray pulsars. 
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The predictions of our analytical model can be tested both observationally and also 
by developing numerical models that allow the relaxation of several of the key assumptions 
made here. For example, within the context of a numerical model one can include the exact 
spatial variation of the shape and strength of the dipole magnetic field, as well as the full 
energy and angle dependence of the electron scattering cross section. Although these effects 
are beyond the scope of the present paper, we nonetheless believe that our main conclusions 
regarding the spectral shape will be verified by more detailed models. 
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Table 1. Input Model Parameters 



Model 


Object 




S 


B{G) 


M (g s-i) 


Te (K) 


ro (m) 


D (kpc) 


1 


Her X-1 


1.45 


1.80 


3.80 X 10^^ 


1.11 X 10^^ 


6.25 X lO''' 


44.0 


5.0 


2 


LMC X-4 


1.15 


1.30 


3.28 X 10^2 


2.00 X IQi** 


5.90 X 10'' 


680.0 


55.0 


3 


Cen X-3 


1.25 


3.71 


2.63 X 10^2 


1.51 X 10^® 


3.40 X 10' 


730.0 


8.0 



Table 2. Computed Model Parameters 



Model 


a 


cr||/(T^ 


cr/a,j. 


J (gs ^ cm ^) 


Tth (K) 


fth/c 




Tth 


'?~max 


''"trap 


1 


0.401 


4.15 X 10-5 


2.93 X 10-* 


1.83 X lO'' 


5.68 X 10' 


0.027 


0.098 


0.067 


1.37 


1.58 


2 


0.319 


4.85 X 10-5 


3.98 X 10-* 


1.38 X 10* 


2.91 X 10^ 


0.026 


0.034 


0.080 


1.12 


1.77 


3 


0.346 


8.30 X 10-5 


4.51 X 10-" 


9.02 X 10^ 


2.48 X 10^ 


0.023 


0.072 


0.066 


1.07 


1.70 



Table 3. Auxiliary Parameters 



Model 


ec (keV) 


CTc (keV) 


dc (keV) 


CK (keV) 


<JK (keV) 


(Ik (s ^ cm 




1 


44.0 


13.0 


25.20 


6.55 


0.30 


4.00 X 10" 


-3 


2 


38.0 


21.2 


38.29 


5.90 


0.22 


1.00 X 10" 


-3 


3 


30.5 


11.0 


21.50 


6.70 


0.27 


1.20 X 10- 


-2 



